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TITULO 

Composicion coherente de senales por correccion focal progresiva. 

SECTOR DE LA TECNICA 
5 Tecnologia fi'sica 

Sistema electronico para Ecografia clmica y Ensayos No Destructivos. 
Otras aplicaciones: radar, sonar y seismologia 

ESTADO DE LA TECNICA 

10 En diversos ambitos de aplicacion, un conjunto de N transductores o array, 
recibe senales de una o multiples fuentes ubicadas a cierta distancia y 
orientacion relativa al array. Esta situacion es comun en radioastronomia, 
radar, sonar, sismologla y, particularmente, en el ambito de la tecnologia 
ultrasonica basada en arrays (ecografia clmica o industrial, en este caso 

15 aplicada a los Ensayos No Destructivos), donde N se encuentra, tipicamente, 
en el rango 32 a 256 o mas. 

En estas situaciones, cada uno de los N elementos del array recibe la senal 
generada por cada fuente con cierto retraso temporal respecto a los restantes 
elementos. La funcion de un conformador de haz es obtener una unica senal a 

20 partir de las N recibidas en la que los retrasos temporales para cada posible 
distancia en una orientacion determinada se compensen, aumentando la 
capacidad de discriminacion del sistema (mayor resolucion, relacion senal/ruido 
y contraste). Repitiendo este proceso en una diversidad de direcciones 
mediante barridos lineales o angulares, se obtiene informacion suficiente para 

25 poder componer una imagen 2D o 3D de las fuentes de senal. En el caso 
particular de los ultrasonidos, las N senales recibidas son ecos producidos por 
discontinuidades en el medio que ha sido insonificado por un emisor, que 
puede ser el mismo array receptor (o un subconjunto del mismo) u otro 
transductor (array o monoelemento). 



En recepcion, es deseable compensar los retardos en todos los puntos o focos 
a lo largo de la Ifnea de propagacibn en una direccion determinada. Esto 
requiere modificar continuamente los retardos en el tiempo, para adaptarlos a 
la llegada de los ecos desde distancias dinamicamente variables a cada 
elemento. Esta tecnica de focalizacion dinamica, proporciona la mayor 
resolucion, contraste y relacion senal/ruido. 

Siendo tecnicamente imposible modificar en forma continua los retardos, es 
esencial que su discretizacion se realice con una excelente resolucion temporal 
para cada elemento del array y para cada distancia, tipicamente del orden de 
1/16 a 1/64 el periodo fundamental de la serial. 

Este hecho ha motivado la busqueda de metodos para la focalizacion dinamica 
en recepcion sencillos y eficientes que, a lo largo de los anos, se han ido 
adaptando a las disponibilidades tecnologicas. Asf, se han descrito tecnicas 
basadas en dispositivos CCD accionados por relojes cuyo periodo determina el 
retardo [Macowski, "Ultrasonic electronic lens with reduced delay range", US. 
Pat. Nr. 4.058.003, Nov. 1977] o en lineas de retardo conmutables [Maslak, US. 
Pat. Nr. 4.550.607, "Phased array acoustic imaging system", Nov. 1985]. Estos 
metodos introducen artefactos de conmutacion, son voluminosos y 
proporcionan una resolucion limitada estando, practicamente, en desuso. 

Otra alternativa, asimismo analogica, se basa en la utilizacion de mezcladores 
y desplazadores de fase que realizan los retrasos temporales en el dominio de 
la frecuencia [Maslak et al., "Dynamic focused linear phased array acoustic 
imaging system", US. Pat. 4.699.009, Oct. 1987]. Uno de los problemas 
asociados a este tipo de tecnicas es la dificultad de su adaptacion a 
transductores de distintas frecuencias. 

Por otra parte, se ha propuesto la utilizacidn de rotadores de fase, basados en 
el algoritmo Cordic [R. Andraka, "A survey of CORDIC algorithms for FPGAs", 
Proc. 1998 ACM/SIGDA VI Int'l Symposium on FPGAs, Monterrey, CA, pp. 191- 
200, Feb. 1998]), y su aplicacion a la focalizaci6n dinamica se encuentra en [O" 
Donnell et al., "Method and apparatus for digital phased array imaging", US Pat. 
Nr. 4.983.970, Ene. 1991] y en [Engeler et al. "Dynamic phase focus for 
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coherent imaging beam formation", US Pat. 5.111.695, May 1992]. El hardware 
involucrado es significative. 

La disponibilidad de tecnologia VLSI permitio nuevas propuestas basadas en 
retardar las senales una vez digitalizadas. Cada elemento del array se asocia a 
5 un canal de adquisicion y procesamiento (focalizacion) independiente. Las 
salidas de los N canales se suman para obtener la serial resultante (suma 
coherente). Para operar en tiempo real, cada canal dispone de dos espacios de 
memoria: uno para datos, donde se almacenan las muestras adquiridas, y otro 
de focalizacion que establece los retardos a aplicar. 

10 En un principio, la resolucion de los retardos se limito al propio periodo de 
muestreo de las senales, determinado por el criterio de Nyquist, alrededor de 
1/3 a 1/5 del periodo de la serial, utilizando registros de desplazamiento de 
longitud variable [Welles II et al., "Architecture for ultrasonic imaging", US Pat. 
4.796.236, Ene 1989]. La baja resolucion temporal produce lobulos de 

15 cuantificacion de retardos que reducen el rango dinamico, contraste y relacion 
senal/ruido de la imagen [D. K. Peterson, G. S. Kino, "Real-time Digital Image 
Resconstruction: A Description of Imaging Hardware and an Analysis of 
Quantization Errors", IEEE Trans. Sonics Ultrasonics, 31, pp. 337-351, Jul. 
1984]. 

20 La opcion de utilizar una mayor frecuencia de muestreo no es aconsejable por 
el mayor coste de los componentes, consumo energetico y volumen de datos a 
almacenar y procesar. Para solventar este problema se han propuesto diversas 
tecnicas que, basicamente, pueden agruparse en dos: 

1) Realizacion de un retardo grueso con una resolucion temporal equivalente al 
25 periodo de muestreo, y otro fino mediante interpolacion de la serial recibida 
entre dos muestras consecutivas, con una resolucion, tipicamente, entre 1/4 y 
1/16 del periodo de muestreo. 

A este grupo pertenecen las patentes: [Odell, "Digital beamforming and filtering 
circuit", U.S. Pat. Nr. 5,268,877, Die. 1993], [D.Lipschutz, "Time multiplexed 
30 digital ultrasound beamformer"US.Pat. 5,469,851, Nov. 1995], [Beaudin et al. 
"Delay generator for phased array ultrasound beamformer", US. Pat. 5.522.391, 
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Jun. 96], [Wright et al. "Method and apparatus for receive beamforming 
system", US Pat. 6.042.547, Mar. 2000] y [M.H. Bae, "Ultrasonic signal focusing 
method and apparatus for ultrasonic imaging system", US Pat. 6.231 .51 1 , May. 
2001]. 

5 Esta tecnica, con diversas variantes para la realizacion ffsica de los retardos 
grueso y fino, incluida la de rotacion de fase mediante algoritmos CORDIC, 
requiere invariablemente una etapa de interpolacion y un selector de muestras 
interpoladas por cada canal, con la consiguiente complejidad y errores 
asociados. 

10 2) Adquisicion directa de las senales en los instantes que corresponden a la 
llegada de la serial de cada foco a cada elemento. Para ello, se controla la fase 
del reloj de muestreo de cada canal con una resolucion temporal equivalente 
entre 1/16 y 1/64 el periodo fundamental de la senal. Este metodo hace 
innecesaria la etapa de interpolacion, eliminando la electronica y las fuentes de 

15 error correspondientes. Las muestras se pueden adquirir sobre una memoria 
FIFO, de la que se extraen para realizar la suma coherente [Chang et 
al., "Bandwidth sampling technique for digital focusing in array imaging 
systems", US. PAT. 5,581,036, Die. 1996]. 

El mayor problema asociado a esta tecnica es la realizacion de los generadores 
20 de reloj de muestreo de fase controlada con la suficiente precision. En 
[O'Donnell et al., "Method and apparatus for fully digital beam formation in a 
phased array coherent imaging system", US Pat. Nr. 4.809.184, Feb. 1989] y 
[Corl, "Digital ultrasound system with dynamic focus", U.S. Pat. 4,974,21 1, Nov. 
1990], se generan m relojes derivados de un reloj maestro de frecuencia m 
25 veces superior a la de muestreo, cada uno de ellos con un desfase 2n/m\ un 
selector del reloj por cada canal, accionado por el contenido de la memoria de 
focalizacion, permite asociar en cada instante uno de los relojes al de 
muestreo. La distribution de multiples relojes y su selection con coherencia 
temporal es un problema, especialmente para las mayores resoluciones 
30 temporales. 

Otra posibilidad es almacenar el reloj de muestreo de cada canal en una tabla 
como una secuencia de 1s y 0s, como se ha sugerido en [J.H. Kim, T. Song, 
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S.B. Park, "Pipelined Sampled-Delay focusing in ultrasound imaging 
systems, Ultrasonic Imaging 9, pp.75-91, 1987], pero no es practica por la 
cantidad de memoria necesaria para almacenar todos los relojes de muestreo 
de todos los canales para multiples lineas de barrido. 

5 Siguiendo otra alternativa, se calculan en tiempo real los instantes de muestreo 
mediante circuitos especializados asociados a cada canal. Esta tecnica se 
describe en [Park et al., "Real time digital reception focusing method and 
apparatus adopting the same", US Pat. 5.669.384, Sep. 1997], [Bae, " Focusing 
delay calculation method for real-time digital focusing and apparatus adopting 
10 the same" US Pat. 5,386.881, Feb. 1998] y [Petrofsky, "Method and apparatus 
for distributed focus control with slope tracking", US Pat. 5.724.972, Mar. 1998]. 
Todas, con diferentes variantes, se basan en el algoritmo del punto medio 
descrito en [Bae et al., An Efficient Real Time Focusing Delay Calculation in 
Ultrasonic Imaging Systems, Ultrasonic Imaging, 16, pp. 231-248, 1994]. 

15 Aunque se trata de una tecnica interesante por sus posibilidades de integracion 
en un ASIC, presenta diversos inconvenientes entre los que destacan: la 
necesidad de programar un gran numero de parametros para cada Ifnea de 
barrido y cada canal, la complejidad del hardware para obtener una buena 
resolucion, menor eficiencia de utilizacion de area de silicio frente al requerido 

20 por sistemas basados en memoria por la baja regularidad del hardware y, 
finalmente, requiere un compromiso entre obtener una alta velocidad de calculo 
o una reduccion del rango de aplicacion. 

El metodo y sistema electronico objeto de la presente invencion pertenece a 
este grupo y resuelve los problemas descritos mediante una tecnica mixta, que 
25 combina una memoria de focalizacion de reducido tamano, en la que solo se 
requiere una fraccion de bit por cada foco y canal, con el calculo en tiempo real 
del instante de muestreo con elevada resolucion mediante circuitos sencillos y 
regulares de facil integracion VLSI. 



30 
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DESCRIPCION DE LA INVENCION 
INTRODUCCION: 

Se describe un conformador de haz que realiza la composicion coherente de 
las senales procedentes de una direccion de propagacion determinada, 
5 recibidas por un conjunto de N transductores receptores o array. La 
metodologia y aparatos descritos son aplicables en campos como imagen 
ultrasonica (ecografia y Ensayos No Destructives), radar, sonar, acustica, 
seismologfa y otros. 

El conformador de haz aquf descrito opera con arrays cuyos elementos pueden 
10 estar distribuidos arbitrariamente. En particular, realiza la focalizacion dinamica 
con deflexion del haz en acimut y elevacion utilizando arrays bidimensionales. 

Basa su operacion en muestrear la serial recibida por cada elemento del array 
en los instantes que corresponden a los de llegada de la senal desde multiples 
focos situados a lo largo de una direccion de propagacion dada. Dentro de los 

15 Ifmites establecidos por los principios que gufan la presente invencion, el 
instante de adquisicion de la serial procedente de cada foco y recibida por cada 
elemento del array se determina con un error absolute respecto al valor 
exacto, inferior a medio periodo de un reloj maestro. La frecuencia del reloj 
maestro puede elegirse para reducir la incidencia negativa de los lobulos de \ \ 

20 cuantificacion de los retardos en la relacion senal/ruido, contraste y rango 
dinamico de la senal resultante. 

La frecuencia de muestreo es independiente de la de la senal, siempre que se 
verifique el criterio de Nyquist, tanto en radiofrecuencia como en banda base, y 
de la frecuencia del reloj maestro. Con tecnicas de sobremuestreo (conversion 
25 AZ), la frecuencia de digitalizacion puede ser la mas baja compatible con. la 
relacion senal/ruido deseada. En particular, puede utilizarse el reloj maestro 
modificado segun el metodo descrito para realizar directamente el muestreo 
AZ. 

Al reducir la frecuencia de muestreo a la mas baja posible, compatible con el 
30 ancho de banda de las senales y/o relacion senal/ruido, se reducen el coste, 
consumo energetico y volumen de datos a procesar. En particular, con el 



8 



metodo que constituye una parte esencial de la presente invention, solo se 
adquieren las muestras necesarias para realizar la suma coherente, sin que 
sea necesario un proceso de interpolation o la adquisicion de information 
redundante. 

5 Ademas, el metodo de conformation que se describe y que constituye otra 
parte esencial de la presente invention, posee unos bajos requisitos de 
memoria de focalizacion, que pueden llegar a ser de una fraction de bit por 
cada muestra adquirida en cada canal. Esta alta eficiencia facilita que el 
metodo pueda aplicarse convenientemente a una diversidad de modos de 
10 operation, como son el barrido lineal o sectorial, con aperturas densas o 
dispersas, y otras aplicaciones en las que se requieren multiples descriptores 
de la focalizacion dinamica para cada configuration de la apertura activa del 
array. 

Dichas caracterfsticas posibilitan la integration conjunta de la memoria de 
15 focalizacion, la de adquisicion de datos, la logica de generation de los relojes 
de muestreo y la de otras funciones auxiliares como apodizacion, apertura 
dinamica y control, de multiples canales, en un unico dispositivo VLSI o 
estandar configurable tipo FPGA, con la consiguiente reduction de coste sin 
perdidas de calidad. 

20 

Descripcion del metodo. 

Por sencillez, la exposition se realiza para un array unidimensional de N 
elementos, no necesariamente equiespaciados, considerando una aplicacion 
de ultrasonidos utilizando la tecnica de pulso-eco en un medio homogeneo con 
25 velocidad de propagation c. El metodo que se describe a continuation se 
denominara genericamente tecnica de Correction Focal Progresiva o CFP, que 
describe la forma en que opera. 

La Figura 1 muestra un diagrama de la geometna general con los 
condicionantes anteriores. El array, de N elementos, numerados 1, 2, 3, /c, 
30 N-1 y A/, se situa sobre el eje de las x, con el origen de coordenadas 

arbitrariamente definido, por comodidad, en el centra del array. Se desea hacer 
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una adquisicion de senal en la direccion angular 0 referente al eje z. Para ello, 
se definen multiples focos, a intervalos AR en la direccion de propagation 
(llegada de la senal). El metodo propuesto combina las senales procedentes de 
cada foco y recibidas por los diferentes elementos del array de forma que se 
5 obtiene una salida que representa la senal que seria proporcionada por un 
transductor del tamano del array focalizado en todos y cada uno de los focos. A 
este proceso se le denomina composition coherente de las senales recibidas y, 
al resultado, senal focalizada dinamicamente con deflexion. 

Los intervalos AR se eligen convenientemente en funcion del periodo de un 
10 reloj maestro T x =1/fx V de la velocidad de propagation c como: 

AR=vcT x /2 (1) 

donde v es un numero entero que expresa el numero de periodos del reloj 
maestro que hay entre dos focos consecutivos en pulso-eco. Por otra parte, los 
intervalos AR deben ser lo suficientemente pequenos como para que se 

15 verifique el criterio de Nyquist: si T s es el intervalo temporal entre muestras, y 
T R es el periodo fundamental de la senal, el criterio de Nyquist se verifica, 
tfpicamente, para T s < Tr/3, considerando el ancho de banda de las senales 
habituales y la adquisicion en radiofrecuencia. Preferentemente, el intervalo 
temporal entre muestras es un multiplo del periodo del reloj maestro, esto es, 

20 T s = sT x , con s entero. 

En general, puede adquirirse un numero m de muestras entre dos focos 
consecutivos, eligiendose preferentemente el valor v en funcion del intervalo 
entre muestras T s - sT x para que v = ms, con m y s enteros. En una 
realization particular en la que m=1, cada muestra coincide con un foco, de 

25 modo que todas las muestras se focalizan dinamicamente. Cuando m>1, 
estrictamente solo se realiza la focalizacion dinamica en los focos, aunque la 
profundidad de foco permite, en general, suponer bien focalizadas las muestras 
intermedias. Una variante del metodo, denominada tecnica de Correction Focal 
Progresiva Variable o CFPV, varia dinamicamente el numero de muestras 

30 adquiridas entre dos focos, reforzando la idea de mantener una buena 
focalizacion dinamica en todas las muestras adquiridas. 
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Por otra parte y, como es bien conocido, el valor [x = T R /T X o relacion entre el 
periodo de la serial recibida y la resolucion temporal del sistema, afecta al 
rango dinamico y contraste de las imageries construidas con multiples haces 
focalizados. En general, se considera aceptable un valor \i > 16 (ver, por 
5 ejemplo, [D. K. Peterson, G. S. Kino, "Real-time Digital Image Resconstruction: 
A Description of Imaging Hardware and an Analysis of Quantization Errors", 
IEEE Trans. Sonics Ultrasonics, 31, pp. 337-351, Jul. 1984]). 

Es importante que la resolucion temporal determinada por T x sea 
independiente del periodo de muestreo T s que se elige en base al criterio de 

10 Nyquist y, por tanto, se relaciona con el periodo de las senates recibidas, T R . 
Cuando el muestreo se realiza en radiofrecuencia se elige, en general, T s * 
Tr/3 o T s « Tr/4 para reducir el flujo de datos a procesar y el consumo 
energetico, esto es, la frecuencia de muestreo mas baja compatible con el 
criterio de Nyquist. En tal caso, con s > 6os> 4, respectivamente, se tiene 

15 que n > 16. Una particularidad de las tecnicas CFP y CFPV es que permiten 
operar a la mayor parte de la electronica a la frecuencia de muestreo, 
reservandose la mas elevada del reloj maestro para muy pocos elementos que 
realizan el ajuste del instante de muestreo. 

Para un foco F, situado en la coordenada polar generica (R, B) y, para el 
20 elemento /c del array situado en la coordenada cartesiana (x k , 0), por el teorema 
del coseno (Figura 1): 



donde L k (R,Q) representa la longitud del camino recorrido por la senal desde el 
foco F situado en (R,Q) al elemento k. El tiempo de vuelo de un pulso 
25 ultrasonico desde su emision arbitrariamente elegido en el origen de 
coordenadas al foco F y de este al elemento k f es: 



En el caso en que el origen de tiempos se situe en una coordenada diferente 
(por ejemplo, cuando el espacio inspeccionado sea iluminado desde un emisor 




(2) 




(3) 
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diferente del propio transductor receptor), el primer termino de la Ecuacion (3) 
debera ser modificado para representar el tiempo de vuelo desde el emisor al 
foco F. 

En general, se denomina T ki al tiempo invertido por la senal correspondiente al 
5 foco F/ en alcanzar el elemento k del array. Esta notacion generica se refiere, 
tanto a la geometna indicada en la Figura 1, como para otras situaciones en las 
que los elementos del array se distribuyen de forma arbitraria, el medio de 
propagacion no es homogeneo o los focos se situan en otras posiciones. No 
obstante, para facilitar la exposicion, se mantiene la geometna de la Figura 1, 
10 el medio homogeneo con velocidad de propagacion c y las posiciones de los 
focos se alinean formando un angulo 0 respecto al array situandose, de 
momenta, a intervalos regulares AR a partir de una distancia inicial R 0 en la que 
se coloca el primer foco F 0 . 

La variation en el instante de llegada al elemento k de senales procedentes de 
15 dos focos consecutivos F h1 y F, situados a distancias R r AR y R, del origen de 
coordenadas, respectivamente, esta dada por la diferencia: 

que, con AR suficientemente pequeno, es: 

dTARfi) 



A7\.CR.,e) = A* 



& i' dR 



(5) 



R = R. 
i 



20 Derivando la Ecuacion (3) respecto a R y sustituyendo el valor de AR dado por 
la Ecuacion (1): 



R. - x, send 

i+- 1 k 

-yjR^ 2R.x k senQ 



~~- (6) 



Esta es una funci6n mon6tona creciente entre R=0 y oo e independiente de 
la velocidad de propagacion, con los siguientes limites: 



12 



lim I 




1- 



senQ 




>0 



(7) 



J 




lim [AT, )=vT 



X 



(8) 



5 



10 

a 



15 



20 



donde, por claridad de notacion, AT k = AT k (R,Q). El valor x k /\x k \ de la Ecuacion 
(7) representa el signo de x k , mientras que 9 interviene con su signo. En 
cualquier caso, el menor valor posible para AT k es cero, independientemente 
de la posicion del elemento k. 

De este modo, los intervalos entre los instantes de muestreo que corresponden 
a dos focos consecutivos estan limitados por los valores expresados en las 
Ecuaciones (7) y (8). En particular, para el caso 9=0 (no se realiza deflexion), la 
funcion crece monotonamente entre vT x /2 y vT x cuando la distancia a la que se 
situan los focos varia de 0 a oo. Cuando se reciben las senales con deflexion 
del haz, el rango de variacion es mayor y, en un caso extremo de deflexion a 
90°, entre 0 y vT x . 

En una posible realizacion, las posiciones de los focos quedan definidas por 
valores comprendidos en el intervalo [0, v], que representan los valores 
incrementales dados por la Ecuacion (6) en periodos de reloj maestro. En este 
sentido, si el numero entero v puede expresarse con b bits, cada retardo 
incremental para cada foco y para cada elemento puede representarse con 
esta misma anchura, independientemente de la posicion ocupada por el 
elemento del array, de la distancia a la que se situe el foco y del angulo de 
deflexion. El error absoluto con el que se determina el instante de muestreo es 
inferior a medio periodo de reloj maestro al redondear al entero mas proximo el 
resultado de AT ki (R,Q)/T x dado por la Ecuacion (6). El inconveniente de esta 
opcion es que requiere un elevado numero de bits para expresar los retardos 
incrementales para cada foco y elemento, con las consiguientes demandas de 
capacidad de memoria. 



Inversamente, si se dispone de b bits para representar cada retardo 
incremental, el rango de valores posibles puede elegirse entre dos situaciones 
extremas: 



(v-2* + l) 



\T x <AT k <vT x 



0<AT k < (2 b - \)T X 



(9a) 



(9b) 



La forma (9a) resulta util para realizar la focalizacion dinamica a partir de una 
distancia minima R min , esto es, para R > R min , mientras que la forma (9b) seria 
preferible para focalizar dinamicamente en las regiones mas proximas al array, 
esto es, para R < R max . Entre uno y otro caso pueden darse rangos intermedios 
mas adaptados a una aplicacion particular, donde R min <R< R max . 

Sin embargo, con la Tecnica de Correccion Focal Progresiva se prefiere la 
forma (9a) que no limita la maxima distancia de validez. Para las regiones mas 
proximas al array en las que no es posible verificar la desigualdad de la 
izquierda de la Ecuacion (9a), la tecnica propuesta proporciona una serie de 
alternativas que seran descritas con detalle mas adelante. Asi, pues, llamando 
en la Ecuacion (9a): 



a = 2 b -I 



puede escribirse en forma mas compacta: 

(v - a)T x < AT k <vT x 



(10) 



(11) 



La Ecuacion (6) siempre verifica la desigualdad de la derecha de la Ecuacion 
(11), de acuerdo con la Ecuacion (8). Para que tambien verifique la de la 
izquierda: 



R - x , send 
k 



11 * 

^R 2 -2Rx h senQ 



vT 



> (v - a)T 



X 



(12) 



Con: 
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2 2Rx send = (r~ x.senQ ) |l + k 



2 2 
x, cos 0 



[r — x.senB p 



la Ecuacion (12) queda como: 



( 



1 + 



1 



1 + 



2 2 
cos 6 

A: 

[ft. - x^senQ 



— >v - a 
2 



(13) 



(14) 



donde, haciendo los cambios de variable: 



se obtiene: 



r = R - x j send 
k 



(15) 
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\ 



i + 



i 



l+ 4"- 2 



v _ 

— >v - a 
2 



Despejando rde la ecuacion anterior se obtiene: 

z (v -2d) 
r > K v>2a 

2^1 av - a 2 



(16) 



(17) 



10 que, deshaciendo los cambios de variable de las Ecuaciones (15) conduce a: 



R>R Q (x k ,e) 



v ^ a = COS0 
>Vav -a 2 



k 



v >2a 



(18) 



donde \ | represents el valor absoluto y considera el caso mas desfavorable de 
signos relatives de x k y 0. La ecuacion (18) representa una cota superior para 
la distancia minima R 0 (x k , Q) a partir de la cual los intervalos temporales para 
muestrear la senal recibida por el elemento k con un angulo de deflexion 0 
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pueden representarse con un numero entero de b bits, con un error inferior a 
medio periodo T x de un reloj maestro. 

Si la Ecuaci6n (18) se verifica para todos los elementos del array a todas las 
distancias y angulos de la region a inspeccionar, el retardo inicial para adquirir 
5 la primera muestra, correspond iente al foco F 0 , en el canal k viene dado por: 



donde R 0 = R 0 (x k , 0) dado por la Ecuacion (18). Obviamente la primera 
adquisicion puede realizarse a una distancia R 1 > R 0 , en cuyo caso el valor de 
To(Xk, 6) se calculara para este nuevo valor. 

10 A partir de este instante, los intervalos de muestreo entre focos presentan un 
error acotado a medio periodo del reloj maestro, pudiendo codificarse para 
cada foco F, de coordenadas (R h 9) en el canal k como: 



donde [ ]n representa la funcion de redondeo al entero mas proximo y AT ki (R„ 
15 e; esta dado por el valor de la Ecuacion (6). Sustituyendo en (20) las 
Ecuaciones (10) y (1 1 ): 



que indica que los intervalos de muestreo pueden codificarse con b bits, 
produciendo un error inferior a medio periodo de reloj maestro. El valor Q ki es el 
20 codigo de correction focal que determina el numero de periodos de reloj 
maestro que debe adelantarse la adquisicion de la serial correspondiente al 
foco / en el canal k. 

La adquisicion de un numero m > 1 de muestras entre focos aumenta la 
eficiencia de memoria al compartir un unico codigo de correccion focal de b bits 
25 entre varias muestras, con una densidad de b/m bits por muestra. Si, por 
ejemplo, b=1 y m=4, se requieren 0'25 bits para focalizar cada muestra 




(19) 




T 

L l x in 



(20) 



0<Q,.<2 b -l 



(21) 
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adquirida. Una buena election para reducir el error en los intervalos de 
muestreo consiste en distribuir uniformemente el adelanto expresado por el 
codigo de correction focal Q ki entre las m muestras. Para esto, en una 
realization preferida, se calculan los valores: 



m 



\<j<m (22) 



donde q= Q ik representa, genericamente, el codigo de correction focal para 
simplificar la notation y [ ]u es la funcion de redondeo. El valor de u(j, q) 
representa el entero mas proximo del numero de periodos del reloj maestro 
entre el primer foco y la muestra j con y=7, 2, m, consiguiendo una 
10 distribution lo mas uniforme posible del adelanto q entre las m muestras. 

El intervalo temporal entre un foco cualquiera con la muestra que le sigue en el 
orden j=m es u(m, q) = v-g, esto es, la muestra j=m se adquiere exactamente 
en la position del siguiente foco, con un adelanto de q ciclos del reloj maestro 
respecto del intervalo nominal v entre focos, de acuerdo con lo expresado por 
15 el codigo de correction focal. 

El intervalo entre muestras, en periodos del reloj maestro, viene dado por: 
Au(l, q) = u(l, q) j = 1, por definicion 
Au(j,q) = u(j\q)-uU-l,q) 2<j<m < 23 > 

En el caso mas habitual en el que el numero de muestras entre focos es mayor 
o igual al maximo adelanto a aplicar entre focos, esto es, q < 2*-1 < m, al 

20 hacer una distribution uniforme de q ciclos entre m muestras, el adelanto g, 
asignado a la muestra j solo puede tomar los valores 0 6 1 respecto al intervalo 
nominal v/m que se da para q=0. Este hecho facilita la realization de la 
distribution uniforme del adelanto q entre un numero de muestras m > q 
mediante una tabla en la que se codifica con 1 bit el valor g(j,q), previamente 

25 calculado segun: 



gU,<l)=v /m-Au(j\q) 1 < j < m y 0<q<2 b -1 



(24) 
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Para clarificar el metodo, considerese un ejemplo en el que m=3, v=12=4 m, 
b-2. Los posibles valores de q son 0, 1, 2 y 3 por estar codificado con 2 bits, y 
v/m=4. La Tabla I muestra los resultados que se obtienen al aplicar, con estos 
valores, las Ecuaciones (22), (23) y (24). Se observa que el adelanto q se 
5 distribuye uniformemente entre las muestras dentro de la resolution de 1 
periodo de reloj maestro para todos los codigos de correction focal. La 
densidad de la codification en este ejemplo es de 2/3 = 0.67 bits/muestra. 

Tabla I. 

10 Ejemplo de distribution uniforme de q adelantos entre m muestras 





v-q 


u(1.q) 


u(2.q) 


u(3,q) 


Au(1,q) 


Au(2,q) 


Au(3,q) 


9(1 .q) 


9(2,q) 


9(3,q) 


q=0 


12 


4 


8 


12 


4 


4 


4 


0 


0 


0 


q=1 


11 


4 


7 


11 


4 


3 


4 


0 


1 


0 


q=2 


10 


3 


7 


10 


3 


4 


3 


1 


0 


1 


q=3 


9 


3 


6 


9 


3 


3 


3 


1 


1 


1 



Para cada valor de q hay m valores posibles g(j,q), cada uno de 1 bit. Por 
consiguiente, todos los adelantos individuates pueden codificarse en una 
memoria de m 2* palabras de 1 bit (en el ejemplo anterior, basta una memoria 
15 de 12x1 bit). 

Para valores R < R 0 (x k , Q) el error absoluto del intervalo entre instantes de 
muestreo puede superar el valor 7>/2. En este caso, las tecnicas CFP y CFPV 
pueden operar de varias formas, destacando las siguientes estrategias: 

1. Activar una funcion de apertura dinamica en la que el elemento k no 
20 interviene hasta que se reciban senales desde una distancia igual o superior 
a Ro(x k , Q). Una de las posibles realizaciones, produce para este elemento 
un numero de muestras nulas M z (k) equivalente a las que se producen a la 
salida del conformador entre el comienzo de la adquisicion T in i y el instante 
T 0 (x k , Q) dado por la Ecuacion (19): 
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M z (k) = 




Jt 



(25) 



donde [ ]t representa la funcion de redondeo por exceso a un numero 
entero. 

Debe destacarse que la relacion D/max(R 0 ), donde D es el tamano de la 
apertura, determina el mmimo valor de F# que puede aplicarse con las 
tecnicas propuestas. En un ejemplo practico se muestra que F# mm =1, 0'6 y 
0'4 para valores de b= 1, 2 y 3, respectivamente, muy competitivos con los 
proporcionados por otras tecnicas. 

2. Anticipar la adquisicion de serial en el elemento k a una distancia R t < 
Rofok, Q) elegida de forma tal que se adquiera una muestra en R 0 (x k , Q). Los 
instantes de muestreo anteriores a este tendran un cierto error superior al 
acotado, pero su efecto en la imagen puede ser limitado. Para ello, el 
instante de adquisicion de la primera muestra dado por la Ecuacion (19) 
queda modificado a: 



donde se expresa que, para el elemento /c, se adquieren M z (k) muestras 
anteriores a la recepcion con error acotado y T S i es el periodo de muestreo 
utilizado, constante en esta region. El valor de T S i puede ser igual al 
periodo de muestreo de salida del conformador 7~s o, preferiblemente, un 
valor inferior para reducir en lo posible los errores en el instante de 
muestreo. 

3. Operar de forma mixta, donde un primer rango de distancias [R 2 , Ri) se 
cubre con aperture dinamica, otro [R 1r R 0 ) con un error limitado como en el 
punto anterior y, finalmente, un tercero [R 0 , <*>) en el que los errores de 
muestreo estan acotados dentro de medio periodo de reloj maestro 
aplicando la correccion focal progresiva. 

En una realizacion preferida que, ademas resulta de aplicacion general cuando 
la geometria del array es diferente, el medio no es homogeneo (por ejemplo, 



T 6 VI (x k #)=T 0 (x k ,e)-M z (k)r, 



SI 



(26) 
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cambios de velocidades de propagacion) o los focos se distribuyen a intervalos 
regulares pero siguiendo una trayectoria diferente a la considerada, se utiliza 
un algoritmo de calculo basado en los siguientes pasos: 

ALGORITMO 1: 

PASO 1 : Se calculan los tiempos de propagacion T ki de la senal 
correspondiente a cada foco F h /= 0, 1, 2, hasta cada elemento k 
del array. El valor de T ki se obtiene en cada caso al considerar la 
geometria del sistema y la(s) velocidad(es) de propagacion. 

PASO 2: Para cada elemento /c, se calculan los intervalos U ki entre instantes 
de llegada de la serial correspondiente a dos focos consecutivos, 
expresados en periodos del reloj maestro, como: 



k f ki 



T -T 
ki kj-l 

T 

l X 



i>\ (27) 

U 



PASO 3: El valor v o intervalo nominal en periodos de reloj maestro entre 
focos se elige de los resultados de la ecuacion anterior como: 

v^maxCC/^.) Vk,i (28) 

PASO 4: Si el numero de bits disponibles para codificar las correcciones 
focales es b, determinar para cada canal k el mdice i=h k a partir del 
cual: 

v- Uk .<2 b -l i>h k (29) 

PASO 5: Para los focos F, con /'< h k , el numero de muestras M z (k), para 
operar con apertura dinamica o adquisicion anticipada en el canal k, 
es: 

M z (k)=(h k -\)m (30) 

PASO 6: Los codigos de correcci6n focal para cada elemento k 
correspondientes al foco / > h k se calculan como: 

Qh=v-U u i>h k (31) 
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que expresan el numero de ciclos de reloj maestro que debe 
adelantarse la adquisicion de la senal correspondiente al foco / en el 
canal k. 

En otra realizacion preferida, durante la adquisicion de serial, la distancia entre 
dos focos consecutivos F h1 y F, se vana en momentos seleccionados 
incrementando el numero de muestras que se adquieren entre ellos. Esta 
variante se denomina Tecnica de Correction Focal Progresiva Variable o 
CFPV, en la que el numero de muestras adquiridas entre dos focos 
consecutivos pasa de ser una constante m a una funcion del numero de foco 
actual, esto es M=m(i). 

La Ecuacion 18 indica que el valor de R 0 crece con el de v. Por su parte, el 
valor mmimo de v esta determinado por la Ecuacion 1 cuando se adquiere una 
muestra por foco, con un interval© temporal entre muestras T s = sT*. Iniciando, 
pues la adquisicion con M=m(1)=1 a partir de una distancia R 0 (1) determinada 
por el valor v = s a la que, para todos los elementos: 



* 0 (1)* 



s -2a 



l4as - a 



= x. COS0 
2 * 



x, send 
k 



(32) 



donde a = 2 b - 1. La adquisicion continua manteniendo M=m(1)=1 hasta 
alcanzar cierto foco a una distancia R=R 0 (2) a partir de la cual se verifica, para 
todos los elementos: 
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25 



* 0 (2)> 



2s -2a 



l^2as-a 2 



■x, cosG 
k 



x, send 
k 



V* 



(33) 



donde comienza un tramo en el que pueden adquirirse 2 muestras por foco, 
esto es, M=m(2)=2, manteniendo el error de muestreo acotado. Cierto tiempo 
despues, se alcanzara otro foco en R-R 0 (3) a partir del cual puede operarse 
con M=m(3)=3 y asf sucesivamente, de forma que, genericamente, a partir de 
un foco situado en R=R 0 (m) a partir del cual pueden adquirirse m muestras, 
verificandose: 
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R Q (m)> 



ms - 2a 



2-/ ams - a 



= x, COS0 
2 k 



+ 



x, send 
k 



V* 



(34) 



En definitiva, con la Tecnica de Correcci6n Focal Progresiva Variable el numero 
de muestras que se adquieren entre focos se incrementa con el tiempo, de 
forma que los codigos de correccion focal puedan seguir expresandose con b 
bits. Debe destacarse que la condicion (34) debe ser verificada por todos los 
elementos para cada valor de m. 

En una realization preferida, se asocia a cada foco F, un campo de 1 bit, 
denominado codigo de alejamiento focal J,. Cuando J ; =0 se mantiene el actual 
valor m(i)=m(i-1) y, cuando Jj=1 se incrementa: m(i)=m(i-1)+1. El algoritmo 2 
proporciona un posible metodo para el calculo de los valores de Jr. 



ALGORITMO 2: 

PASO 1: En funcion de la distancia minima establecida por la aplicacion, R mm , 
determinar el mayor valor m=m 0 para el cual se verifica la Ecuacion 
15 (34). 

PASO 2: Inicializacion: m=m 0 . Variables auxiliares: f=1 (numero actual de 
foco), i=1 (numero actual de muestra). 

PASO 3: Repetir: 

a) Calcular R 0 (m) dado por la Ecuacion (34) y el valor AR = cTs/2, 
20 siendo c= velocidad de propagacion, Ts = periodo de 

muestreo. 

b) Si R mi „ + (i+m) AR > (m+1)R 0 

{ J f = 1; m=m+1;} 
en otro caso { J f = 0;} 
25 c) f=f+1; 

d) i=i+m; 

hasta que / > Ns, siendo N s = numero de muestras a adquirir. 
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El numero de codigos de alejamiento focal iguala al de focos para cualquier 
adquisicion, por lo que puede realizarse facilmente asociando a cada codigo de 
correccion focal Q w el correspondiente codigo de alejamiento focal J ki = Jj en 
5 cada canal k. 

Este esquema resulta adecuado ya que, en general, las senales 
correspondientes a cada foco llegan a los elementos en momentos diferentes y 
su realization es sencilla. De este modo, la memoria de correccion focal tiene 
una anchura de palabra de b+1 bits, repartidos en dos campos: b bits para el 
1 0 codigo Q k} y 1 bit para el J,. 

Una alternativa, que proporciona una mayor eficiencia, almacena en una 
memoria separada los codigos de alejamiento focal J h cuyo contenido es 
compartido por un conjunto de canales. Para resolver el problema de no 
simultaneidad en la adquisicion de los focos, cada canal dispone una pequena 
15 cola donde se introducen los codigos J, obtenidos automaticamente de la 
memoria que los almacena. Cada canal utiliza el valor de J, en cabeza de la 
* cola cada vez que alcanza un nuevo foco. Asl, una unica memoria central 

almacena un codigo J, por foco y unos pequenos FIFOs distribuidos facilitan la 
asignacion de estos codigos comunes a cada canal cuando lo necesita, como 
20 se expone con mas detalle en la description del sistema. 

Por otra parte, puesto que la distancia entre focos puede modificarse 
dinamicamente, en la mayoria de las aplicaciones practicas basta una 
codification de la correccion focal con b-1 bit, lo que proporciona una alta 
eficiencia en el uso de memoria y una simplification de los circuitos de 
25 correccion. En particular, es inmediato realizar la distribution uniforme de 0 6 1 
periodo de reloj maestro entre m muestras. 

Con lo expuesto, la Tecnica de Correccion Focal Progresiva, con su opcion 
Variable, proporciona un metodo que realiza la conformacion dinamica de 
haces en recepcion con una elevada resolucion temporal dentro de un gran 
30 rango de distancias. Adicionalmente, carece de redundancias en la adquisicion, 
no requiere procesos de interpolacion de muestras, permite mantener la tasa 
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de muestreo mas baja compatible con el criterio de Nyquist y, por consiguiente, 
tiene un bajo consumo energetico. Su realizacion es compacta, debido tanto a 
un balance adecuado entre el calculo en tiempo real de las correcciones y el 
contenido de unas memorias precalculado en tiempo no real, como a que la 
5 codificacion de esta information es muy eficiente. 

Todo esto permite configurar un sistema versatil de conformacion de haces en 
recepcion con focalizacion dinamica, deflexion, apertura dinamica y 
apodizacion, cuyos detalles constructivos se exponen a continuation. 

Descripcion del sistema. 

10 El sistema objeto de la presente invention se basa en la tecnica de Correction 
Focal Progresiva descrita anteriormente, incluyendo o no la option Variable, y 
realiza la focalizacion dinamica con deflexion de las senates recibidas por un 
conjunto de transductores o array, incluyendo funciones de apodizacion y de 
apertura dinamica. El sistema se describira con suficiente detalle para que 

15 pueda ser entendido y reproducido por personal con conocimientos de 
electronica digital. 

Con referencia a la Figura 2, un sistema de conformacion de haces en 
recepcion consta, en el caso mas general, de un array (10) compuesto de N 
transductores elementales (11) o, simplemente, elementos, que reciben 
20 senales del exterior, numerados 1, 2, .../c, N. Frecuentemente, estos 
elementos son, a su vez, emisores de serial para conformar un haz en emision, 
pasando a ser receptores una vez que esta finaliza. En la figura se han omitido 
los dispositivos necesarios para conformar el haz en emision. 

Las senales recibidas por los diferentes elementos son acondicionadas y 
25 amplificadas por los dispositivos AAS (12), que entregan las N senales 
analogicas a 1f a 2 , a k , a N a otros tantos conversores analogico-digitales 
A/D (13). Los dispositivos AAS suelen incorporar, ademas de la funcion de 
amplificacion, otras que incluyen la compensacion de la atenuacion con la 
distancia mediante un control ganancia-tiempo, el filtrado y, en ocasiones, el 
30 desdoblamiento en las componentes en fase y cuadratura. En este ultimo caso, 
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cada senal analogica a k esta formado por el par (i k , q k ) que representan, 
respectivamente, a las componentes en fase y en cuadratura. 

Una vez acondicionadas, las senales analogicas son digitalizadas por los 
conversores A/D. En el caso de senales en banda base con pares (/*, q k ) 
5 pueden utilizarse dos conversores A/D por cada senal o multiplexarse un unico 
conversor A/D entre dos senales, las muestras pares para la componente i k y 
las impares para la q k o viceversa. 

Otra alternativa es la utilizacion de conversores AS, con los que las senales 
analogicas a k se representan mediante una secuencia continua de bits. En este 
10 caso la tecnica de Correccion Focal Progresiva y sus variantes resulta 
especialmente util al ajustar los instantes de muestreo a los de llegada de la 
senal a cada elemento con la resolucion de un reloj maestro cuyo periodo se 
adapta, ciclo a ciclo, al requerido por tales conversores AE para obtener una 
determinada relacion senal/ruido. 

15 Las tecnicas mencionadas son conocidas por los expertos en el tema y de 
aplicacion general. Las diferencias aparecen principalmente en el dispositivo 
CONF (14) o conformador de haz, cuya funcion es la de combinar las N 
secuencias numericas e 1t e 2 , e kf e N proporcionadas por los conversores 
A/D para crear una unica secuencia de salida r que representa, en forma 

20 digital, a la senal recibida por el array focalizada, idealmente, en todos los 
puntos de una direccion determinada. 

Cada muestra entregada por los conversores A/D es una palabra de w bits 
siendo, tfpicamente, 6 < w < 12 para el caso de conversores de muestreo 
instantaneo y 1 < w < 2 para los de tipo AZ, aunque el sistema puede operar 
25 con otras anchuras de palabra. 

Es, precisamente, la forma en que el dispositivo CONF realiza la funcion de 
combinar coherentemente las N entradas e k para generar la salida r lo que 
constituye el aspecto diferencial de la presente invencion. 

El dispositivo CONF (14) que es objeto de la presente patente, basa su 
30 funcionamiento en la tecnica de Correccion Focal Progresiva que constituye 
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uno de sus aspectos claves, aportando diversas ventajas frente a otros tipos de 
realizaciones. 

Asf, segun otros metodos resenados en la description que se hizo del estado 
del arte, cada muestra de la secuencia numerica e* suministrada por k-esimo 
conversor A/D es sometida a un proceso de retardo como indica la Figura 3. 
Cada elemento se asocia a un canal de procesamiento individual (21). El 
contenido de una memoria de focalizacion dinamica (25) es utilizado por un 
dispositivo de control (24) para generar senales que producen un retardo 
grueso (22) y otro fino (23). El retardo grueso tiene, tfpicamente, la resolution 
de un periodo del reloj de muestreo y el retardo fino una fraction del mismo. 
Estos retardos se obtienen con diversos metodos que han sido descritos en la 
literatura especializada y son bien conocidos. 

El contenido de la memoria de focalizacion dinamica (25) es tal que, para cada 
muestra y para cada canal, los retardos introducidos alinean en el tiempo las 
salidas f 1f f 2 , f N que corresponden a los focos sucesivos F 1f F 2 , F N . Estas 
senales son sumadas en (20) para obtener la secuencia r focalizada 
dinamicamente. 

La mayor ventaja de estas tecnicas reside en que solo es necesario un reloj de 
muestreo, comun a todos los canales. Pero adolece de diversos inconvenientes 
derivados de la necesidad de realizar procesos de interpolation para obtener 
una suficiente resolution temporal en los retardos introducidos a las senales, 
con el consiguiente aumento de complejidad electronica y la posibilidad de 
introducir errores. Ademas, requieren una considerable cantidad de memoria 
para representar los retardos grueso y fino. 

El sistema objeto de la presente invention opera de una forma diferente, 
muestreando directamente las senales en el instante en que Megan a cada uno 
de los elementos del array, procedentes de cada uno de los focos situados a lo 
largo de una direction determinada. Al evitar la introduction de retardos a las 
senales, se hacen innecesarios los procesos de interpolation, eludiendo 
posibles errores y eliminando los circuitos electronicos correspondientes. 
Asimismo, solo se adquieren las muestras necesarias para realizar la 
composition coherente, eliminando redundancias. 
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En principio, el coste del metodo objeto de esta patente es que cada conversor 
A/D debe operar con un reloj de muestreo diferente, de frecuencia no uniforme. 
Como se demuestra con la descripcion que sigue, la realizacion de esta 
caractenstica es sencilla a partir de la metodologfa enunciada como tecnica de 
5 Correccion Focal Progresiva, con la opcion Variable, pudiendo integrarse 
multiples canales en un unico dispositivo tipo FPGA estandar y comercialmente 
disponible. 

Con referencia de nuevo a la Figura 2, el metodo propuesto requiere que el 
dispositivo de conformacion de haz CONF (14), proporcione un reloj de 

10 muestreo independiente y de frecuencia no uniforme CK 1t CK 2 , CK k , CK N 
a cada conversor A/D (13), de forma que se adquieran las senales analogicas 
a 1t a 2 , a ks a N proporcionadas por los acondicionadores de senal AAS 
(12) en los instantes que corresponden a los de su llegada desde cada foco a 
cada elemento (11) del array (10). Las muestras obtenidas por los N 

15 conversores A/D (13) son procesadas para obtener su suma coherente, es 
decir, un valor de salida r que es la suma de las N muestras correspondientes a 
la senal recibida de un mismo foco por los N elementos del array. 

En la presente invencion, esta metodologfa se aplica en forma modular, 
asociando un canal de procesamiento a cada elemento del array. La Figura 4 
20 muestra la estructura de uno de los canales (30) que procesa la senal recibida 
por el elemento k (11) del array. La estructura de todos los canales es identica 
a la mostrada. 

Como en otros casos, la senal es amplificada y acondicionada por el dispositivo 
AAS (12) por metodos conocidos, que incluyen dispositivos electronicos para 

25 amplificar, filtrar y manipular de forma analogica la senal. A la salida del AAS se 
obtiene la senal a k que es digitalizada por el conversor A/D (13). Las muestras 
son adquiridas en los instantes definidos por un generador de reloj (70) que 
produce una secuencia de pulsos CK k no uniformemente espaciados, cuya 
descripcion detallada se realiza posteriormente. El generador de reloj es 

30 habilitado por la senal HFD generada en el dispositivo HFOC (80) que marca 
los instantes en los que el canal comienza y termina la adquisicion de senal con 
focalizacion dinamica. 
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La funcion del generador de reloj es la de garantizar que los instantes de 
muestreo determinados por CK k corresponden a los de la llegada de la senal 
desde cada uno de los focos al elemento /c, y opera segun los principios 
descritos anteriormente como Tecnica de Correccion Focal Progresiva, con su 
posible opcion Variable. 

Con esto, la secuencia de muestras e* obtenida a la salida del conversor A/D 
corresponde, precisamente, a los valores ordenados de la senal recibida por el 
elemento k procedente de cada uno de los focos. Estos valores son 
procesados por un bloque de apodizacion y apertura dinamica APD (60), que 
permiten mejorar las caracteristicas del haz y que produce la secuencia f k . Al 
existir una correspondencia 1:1 de las muestras a las f k , esta secuencia 
tambien esta ordenada con referenda a los instantes de llegada de la senal al 
elemento k desde cada foco. Los sucesivos valores de f k se almacenan 
transitoriamente en una cola o memoria FIFO que absorbe la diferencia en los 
instantes en que se producen muestras correspondientes a un mismo foco en 
diferentes canales. 

Elemento clave es, pues, el generador de reloj de muestreo (70), que produce 
la senal de digitalizacion CK k del k-esimo conversor A/D, de forma que las 
muestras obtenidas correspondan a los valores de la senal recibida por el 
elemento k desde cada uno de los focos. La Figura 5 muestra el diagrama de 
principio de una de las realizaciones preferidas del generador del reloj CK k , 
cuyo funcionamiento se basa en el metodo descrito anteriormente como 
Tecnica de Correccion Focal Progresiva, para el caso en que las correcciones 
focales se codifican con 6=1 bit. 

El circuito mostrado en la Figura 5 contiene solo los elementos basicos 
omitiendose otros dispositivos auxiliares, que son necesarios para el correcto 
funcionamiento, en favor de una mayor claridad de exposicion. Tiene un 
registro REG-A (71a) en el que se carga un valor relacionado con el parametro 
v o intervalo nominal entre focos en periodos del reloj maestro. Su salida V 
consta de multiples Ifneas que determinan el numero de etapas que componen 
el registro de desplazamiento de longitud programable SHR-A (72). Este 
dispositivo obtiene una salida p que es una version retrasada una cantidad de 
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V periodos de reloj c de su entrada CK k . Este y otros dispositivos en el 
esquema mostrado operan sincronamente con el reloj maestro cde periodo Tx- 
La serial de entrada HFD habilita la operacion del registro de desplazamiento 
(72) y del flip-flop (77), de forma que estos dispositivos mantienen su estado 
5 actual independientemente del valor a su entrada cuando HFD-O. Cuando 
HFD=1, se habilita la focalizacion dinamica, interrumpiendose cuando HFD 
retorna a cero. Esta senal se genera en el dispositivo HFOC (80), Fig. 4. 

La senal p es, a su vez, retrasada 1 periodo de reloj maestro mediante el flip- 
flop FFA (73), produciendo la senal pr. El multiplexor MUX (74) proporciona una 
10 salida s que es, bien la serial p producida por el registro de desplazamiento o 
su version retrasada pr, en funcion del estado de la linea de seleccion Q: 

s=psiQ=0; s=prsi Q=1 (35) 

El estado del selector Q lo determina el contenido de una memoria MEM (75) 
en la direccion F proporcionada por el contador CNT (76). En el caso actual, Q 
15 es una simple linea (1 bit). 

La memoria se carga previamente, mediante circuitos no mostrados, con los 
codigos de correccion focal Q ki , donde k es el numero del canal actual y / es el 
numero de orden del foco, con / > 1. Los circuitos para la adquisicion de senal 
en el primer foco (/=0), que se describen posteriormente, determinan la 
20 activacion de la senal HFD. 

La salida del multiplexor (74) se registra en el flip-flop FFB (77), cuya salida 
CK k es el reloj de muestreo del conversor A/D del canal actual k (13), Figura 4, 
acciona al contador CNT (76) y es la entrada del registro de desplazamiento 
SHR-A (72). 

25 Para su funcionamiento, el registro de desplazamiento SHR-A (72) se inicializa 
a 'O 1 , el flip-flop FFA (73) a , 0 , J el flip-flop FFB (77) a T y el contador (76) con la 
direccion de memoria donde se almacena el codigo de la correccion focal 
correspondiente al foco F 1 (i=1). Los circuitos auxiliares para establecer estos 
estados iniciales se omiten en la figura. Asi, los valores iniciales de las 

30 diferentes senales son: CK k =1, p-0 t pp=0 y s=0 (independientemente del 
estado de Q). 
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A partir de este estado inicial, en el ciclo de reloj maestro c que sigue a la 
activacion de la senal HFD=1, se pone CK k =0 al encontrarse en este estado la 
entrada s del flip-flop FFB (77) y el reloj habilitado por su entrada ce. Durante 
los V primeros periodos del reloj c se mantienen las senales p=0 y pr=0. En el 
ciclo siguiente se pone p=1, reflejando que el valor inicial CK k =1 se ha 
propagado por el registro de desplazamiento (72), pero todavfa se mantiene 
pr=0. En este mismo ciclo, si Q=0, es s=pr=0, pero si Q=1, es s=p=1. Por 
consiguiente, en el siguiente ciclo, la salida CK k =1 solo si Q=1; tambien en 
este ciclo se pone pr=1 por lo que, si Q=0 se pondra la salida CK k =1 en el ciclo 
siguiente. Una vez que se hace CK k =1 el proceso se repite. 

Por consiguiente, el numero de ciclos n que transcurren desde el anterior valor 
CK k =1 y el proximo lo determina el valor Q de la siguiente forma: 

Si Q=0: n=V+1 

Si Q=1: n=V+2 (36) 

Acorde con la Ecuacion (20), cuando Q=0 el numero de ciclos entre focos es n 
= v y, cuando Q=1, debe ser n= v-1. Con esto: 

V=v-2 (37) 

que indica el valor inicial a cargar en el registro REG-A (71a) en funcion de la 
distancia nominal entre focos v expresada en periodos del reloj maestro. 

Cada vez que se produce un reloj de muestreo (CK k =1), se incrementa la 
direccion F proporcionada por el contador CNT (76), con lo que la memoria 
MEM (75) coloca en su salida el codigo correspond iente a la correccion focal a 
aplicar en el siguiente foco. 

Se observa que el conjunto FFA (73)-MUX (74)-FFB (77) configuran un registro 
de desplazamiento de longitud variable segun el codigo Q, que puede utilizarse 
como base para las realizaciones descritas a continuacion, aunque hay 
alternativas mas eficientes. Asi, la Figura 6 muestra un esquema de otra 
realizacion preferida que no limita el ndmero de bits con el que se representan 
los codigos de correccion focal Q. En este caso el conjunto FFA-MUX-FFB 
anterior se ha sustituido por un registro de desplazamiento SHRB (78) cuya 
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longitud programable la determina el complemento del codigo de correccion 
focal Q obtenido mediante los inversores (75a) que producen las salidas A/Q. 

En este caso el intervalo entre los instantes de muestreo en dos focos 
consecutivos en el canal actual se determina por el valor multibit Q entregado 
5 por la memoria (75). Aqm el registro de desplazamiento SHRB (78) se inicializa 
a T. 

El funcionamiento es, por lo demas, analogo al caso anterior con la salvedad 
de que el numero de ciclos esta determinado por los valores V y A/Q segun: 

n = V + NQ (38) 

10 Teniendo en cuenta que NQ=2 b -1-Q, el valor inicial a cargar en el registro REG- 
A (71a) en funcion del intervalo nominal entre focos en periodos del reloj 
maestro es: 

V=v-2 b +l (39) 

En una realizacion practica, los inversores (75a) son innecesarios si en la 
15 memoria de codigos de correccion focal (75), se almacenan los valores 
complementarios A/Q en lugar de los de Q, manteniendose en la descripcion 
por coherencia en la notacion. Asimismo, en las realizaciones practicas de los 
circuitos mostrados en las Figuras 5 y 6 se incluye, ademas, logica auxiliar que 
evita la situacion en la que todas las salidas estan en estado '0* de la que, 
20 evidentemente, no saldrian los circuitos. 

La logica descrita permite focalizar cada una de las muestras adquiridas. 
Cuando se prefiera realizar una correccion focal periodica, compartiendo un 
mismo codigo de correccion focal Q cada m muestras, los circuitos anteriores 
se modifican como indica la Figura 7. El valor de m se carga en un registro 

25 RCM (71b) y el valor almacenado M es utilizado por una logica de control 
CTRL (90) junto con el valor Q del codigo de correccion focal entregado por la 
memoria (75) para generar las senales ng y cef. La primera controla la longitud 
del registro de desplazamiento SHRB (78) haciendola igual a ng+1, mientras 
que cef habilita al contador CNT (76) que proporciona la direccion F a la 

30 memoria de correccion focal (75). 
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El registro de desplazamiento SHRA (72) es de longitud programable en 
funci6n del valor V programado en el registro REG-A (71a), con un numero de 
etapas igual a V+1. El registro de desplazamiento SHRB (78) es tambien de 
longitud programable con un numero de etapas igual a ng+1. 

El numero de periodos de reloj maestro que median entre dos ciclos 
consecutivos CK k es igual al numero de etapas interpuestas entre la entrada de 
SHRA (72) y la salida de SHRB (78) que, acorde con lo anterior, es: 

n = V + ng + 2 (40) 

Programando en el registro REG-A (71a) un valor 

V =u/m-3 (41) 

donde v/m representa el intervalo nominal entre muestras en periodos del reloj 
maestro, la Ecuacion (40) queda como: 

n —v I m + ng — I (42) 

La logica de control CTRL (90), mostrada en detalle en la Figura 7, controla el 
intervalo n mediante la serial ng, distribuyendo uniformemente el adelanto de 
periodos de reloj maestro Q entre las m muestras que se adquieren entre dos 
focos consecutivos. 

Para ello, el contenido CM de un contador modulo m, CMOD (91), que se 
actualiza con cada reloj CK k por su entrada de habilitacion ce, proporciona 
parte de la direccion DM de la memoria MDM (92). La otra parte de la direccion 
es el codigo de correction focal Q. La direccion proporcionada a la memoria 
MDM es: 

DM(J,QHJ-l>2 b +Q (43) 

donde j= CM+1 indica el numero de orden de la muestra actual. Programando, 
en tiempo no real, el contenido de las direcciones DM(j, Q) de la memoria MDM 
con los valores g(j,Q) que resultan de aplicar las Ecuaciones (22) a (24), la 
salida g representa el adelanto individual a aplicar a cada una de las muestras 
durante la adquisici6n. El inversor (93) obtiene la serial ng, complementaria de 
la g, de modo que, conforme a la Ecuacion (42): 
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Si g=0 => ng-1 => n = v/m 
S/ g=7 => r?g=0 => n = v/m - 1 (44) 

con lo que se produce un adelanto de 1 ciclo de reloj cuando g(j,Q) =7 y se 
mantiene el intervalo nominal v/m cuando g(j,Q) = 0, tal y como quedo descrito 
5 en el metodo de Correccion Focal Progresiva. 

Por su parte, cada vez que el contador CMOD (91) rebosa (esto es, conmuta 
del valor m-7 a 0), se produce un pulso cef que habilita el incremento del 
contador de focos CNT (76), proporcionando una nueva direccion F+1 a la 
memoria MEM (75). Asf, el proceso se repite con un nuevo valor de Q 
10 suministrado por la memoria MEM (75) con el contenido de la nueva direccion. 

Deben hacerse algunas observaciones adicionales. La primera es que, 
evidentemente, los inversores (93) se eliminan si en la memoria se programan 
los valores complementarios de los g(j,Q) dados por la Ecuacion (24). Otra es 
que la memoria MDM (92) es de pequeno tamano; por ejemplo, una memoria 
1 5 de 256x1 bits puede servir para valores m< 64 con b=2, o para m< 32 con b=3. 

A partir de este esquema se deriva la variante denominada como tecnica de 
Correccion Focal Progresiva Variable. Para ello, el registro RCM (71b) se 
sustituye por un contador, de forma que el numero de muestras m entre focos 
puede irse incrementando en el tiempo a partir de un valor minimo inicial m? 
20 que, tlpicamente, puede ser igual a 1 . 

Tal y como fue descrito, segun la Ecuacion (18), la distancia minima R 0 a partir 
de la cual el error absoluto del instante de muestreo se mantiene acotado en 
±T x /2 crece con el valor del intervalo entre focos v, expresado en periodos T x 
del reloj maestro. Al reducir la distancia entre focos en la region mas proxima al 

25 array mediante la programacion de un valor bajo para mi, disminuye el valor de 
R 0 . Los instantes en que se incrementa el contador RCM (71b) vienen 
determinados por los instantes en que las senales proceden de focos situados 
a suficiente distancia como para que pueda adquirirse una muestra mas entre 
focos en todos los canales. El numero de orden de los focos en que esto ocurre 

30 se calcula en tiempo no real y se codifican en una tabla comun a todos los 
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canales, a un subconjunto de canales o a un canal individual, estableciendo un 
balance entre eficiencia de uso de memoria y modularidad. 

En una realizacion preferida y, en tiempo de adquisicion, los focos en los que 
se incrementa el contador RCM (71 b) se determinan por el valor de un codigo 
de alejamiento focal J, expresado con 1 bit, con / = 1, 2, ... asociado a cada 
foco: el valor actual de m contenido en el contador RCM (71b) se incrementa si 
Jj = 1, manteniendo su valor actual si J, = 0 hasta el siguiente foco, donde se 
repite este proceso. El incremento de RCM s6lo se habilita, por consiguiente, 
cuando se produce el reloj de muestreo para adquirir la serial correspondiente 
a un foco. 

En una realizacion particular, la memoria de correccion focal (75) proporciona 
dos campos: uno, J ki = J,- de un bit, indica si se incrementa el contador RCM 
(71b) en este foco y, el otro Q w de b bits, el adelanto a distribuir entre las m 
muestras actuales. En la figura se muestra la senal J=J-, con una Ifnea de 
trazos. En este caso, preferentemente, se utilizan codigos de correccion focal 
Q ki con b=1 bit, puesto que puede controlarse de forma separada la distancia 
entre focos controlando el numero de muestras m intercaladas. En tal caso, la 
distribucion de los adelantos expresados por Q ki son de 0 6 1 ciclo de reloj 
maestro, realizandose en la muestra cuyo orden es la mitad del contenido de 
RCM (71b), aunque son posibles otras alternativas mas simples con poca 
incidencia en la calidad de la imagen. 

Sin embargo, para mejorar la eficiencia, en otra realizacion alternativa que 
tambien se muestra en la Figura 7, los codigos de alejamiento focal J, 
proceden de una memoria separada (MJ (59), Figura 10) cuya salida es 
compartida por multiples canales por la entrada J a una cola FIFO local (79). 
En tiempo de adquisicion, se realiza una lectura de un nuevo codigo de 
alejamiento focal J, en la memoria global MJ (59) de la Figura 10, cuando 
alguno de los canales a los que sirve detecta que su FIFO local (79) esta vacio. 
El c6digo J-, suministrado se escribe simultaneamente en todos los FIFOs 
locales (79) para ser utilizados por cada canal cuando adquieren senal de un 
nuevo foco, esto es, se activa su senal cef. En este momento, el valor en 
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cabeza de la cola Jc se utiliza para incrementar o no el valor de m contenido en 
el contador RCM (71b). 

Debe destacarse que las colas FIFO locales (79) son de pequeno tamano, 
tipicamente inferiores a 16x1 bit porque solo ha de almacenar un numero de 
codigos de 1 bit equivalente a la maxima diferencia de focos adquiridos en 
diferentes canales dentro de un mismo modulo. De este modo, con esta 
alternativa se consumen menos recursos que con la anterior en la que los 
codigos de alejamiento focal J, se almacenaban en la memoria de correccion 
focal (75) de cada uno de los canales. 

Con la Tecnica de Correccion Focal Progresiva y Variable la utilizacion de 
codigos de correccion focal se hace mas eficiente, distribuyendose entre un 
numero bajo de muestras al principio de la adquisicion (en el Ifmite, 1 bit por 
muestra) y entre un numero mayor de muestras conforme la distancia de los 
focos al array aumenta (una fraccion decreciente de un bit por muestra). 



Por otra parte, es preciso determinar con precision el instante en el que se 
adquieren las muestras que corresponden al primer foco, de coordenadas (R 0 , 
Q), a partir del cual opera la correccion focal segun la metodologfa descrita. 
Este instante esta determinado por la Ecuacion (19) para la geometria de la 
Figura 1. Debe hacerse la observacion que el foco inicial puede situarse a una 
distancia mayor en funcion de las necesidades de la aplicacion. Tambien se 
tendra en cuenta la utilizacion o no del un numero variable de muestras entre 
focos para determinar este valor inicial. 

Esta funcion se realiza en el bloque HFOC (80) de la Figura 4, expuesto con 
mas detalle en la Figura 8. Una realizacion utiliza, simplemente, un contador 
accionado por el reloj maestro que se carga inicialmente con el valor N A (k): 



T 0 (k) 

L T x . 



(45) 

n 



donde T 0 (k) representa el intervalo entre el origen de tiempos hasta que la 
senal correspondiente al foco F 0 llega al elemento k, determinado por la 
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geometrfa y las velocidades de propagaci6n, o como expresa la Ecuacion (19) 
para la configuration de la Figura 1. 

Las senales que corresponden al primer foco F 0 se adquieren, en el canal k, 
cuando rebosa este contador, habilitando, a partir de este momento, la 
actuation de los circuitos de correccion focal progresiva. El error absoluto del 
instante de adquisici6n de este primer foco es inferior a medio periodo de reloj 
maestro en todos los canales. Sin embargo, para mantener esta resolution 
temporal, este contador debe ser accionado a la frecuencia del reloj maestro, 
normalmente elevada y muy superior a la de muestreo. 

Para evitar este inconveniente, una realizacidn preferida, descompone el 
contador de ciclos N A en dos partes, como se muestra en la Figura 8, de forma 
que (omitiendo el indice de canal k): 

N A = K N A1 +N A2 + 1 (46) 

con K, N A i y N A2 enteros. El valor N A1 se carga en el contador CNT-A (82) 
mediante la serial Ida. Cuando se habilita su funcionamiento mediante la 
activation de la serial cea, se acciona el contador con un reloj cb cuya 
frecuencia es K veces inferior a la del reloj maestro, c, evitando las dificultades 
mencionadas. 

Tras N A1 ciclos cb, equivalentes a K N A1 periodos T x , el contador alcanza el 
valor terminal que alza la serial tca=1. En el siguiente ciclo de reloj cb se activa 
el flip-flop FF tipo RS (83), colocando la serial aj=1. Este valor se propaga por 
el registro de desplazamiento SHRG (84) cuya longitud la establece el valor 
programado para N A2 , accionado por el reloj maestro c. Transcurridos N A2 + 1 
periodos de reloj maestro, el valor aj=1 alcanza la salida del registro de 
desplazamiento, con lo que la salida pasa a HFD=1, habilitando a los circuitos 
de focalizacion dinamica una vez transcurridos N A periodos de reloj maestro, 
como indica la Ecuacion (46). De este modo, solo el registro de desplazamiento 
SHRG (84) debe operar a la elevada frecuencia del reloj maestro. 

Tal estado permanece estable hasta que se haya realizado un numero N F de 
correcciones focales, programado en el contador CNT-F (86). Cada correccion 
focal es indicada por la serial de entrada cef, que proporcion la logica de control 
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CTRL (90) de la Figura 7 cuando se toman mutiples muestras entre focos, o es 
el reloj de muestreo CK k si solo se adquiere senal en los focos (Figuras 5 y 6). 

Cuando CNT-F alcanza la cuenta terminal, se alza la senal tcf=1, provocando 
la puesta a cero del flip-flop FF (84) por su entrada R. Al realizarse la funcion 
5 AND (81), la senal HFD retorna a 0, finalizando el proceso de correccion focal 
progresiva. 

No obstante, la adquisicion de senales puede continuar con periodo de 
muestreo constante e igual al nominal v/m hasta completar el numero de 
muestras a adquirir, A/ Sl programado en el contador CNT-S (88). El 
10 rebosamiento de este valor activa la senal FIN que termina por completo la 
adquisicion de senal en este canal. 

Es importante observar que, finalizado el proceso de correccion focal 
progresiva al hacerse HFD-0, la direccion F de la memoria de codigos de 
correccion focal (75) queda congelada en el ultimo valor alcanzado por el 
15 contador (76) de las figuras 5, 6 y 7, lo que posibilita iniciar automaticamente la 
siguiente adquisicion con los valores de correccion focal almacenados a 
continuacion, haciendo innecesario programar un nuevo valor inicial de F. De 
este modo es posible operar con diversas estrategias, como por ejemplo: 

a) Generacion de imagen sectorial, donde cada adquisicion se realiza con un 
20 angulo de deflexion diferente. Los codigos de correccion focal que 

corresponden a los angulos de deflexion sucesivos se almacenan 
secuencialmente en memoria, evitando tener que programar la direccion inicial 
F de la memoria de codigos de correccion focal antes de cada adquisicion. A 
diferencia de otras tecnicas, esta permite mantener el error de muestreo 
25 acotado y constante para todos los angulos de deflexion y a todas las 
distancias. 

b) Realizacion de barrido lineal u otras formas de adquisicion donde una 
subapertura activa compuesta por un subconjunto de los elementos del array, 
cambia en cada adquisicion, eliminando ciertos elementos e introduciendo 

30 otros. Analogamente, los codigos de correccion focal que corresponden a las 
configuraciones sucesivas se almacenan secuencialmente en memoria, 
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evitando tener que cargar la direccion inicial F de la memoria de codigos de 
correccion focal antes de cada adquisicion. 

c) Estrategias mixtas, en las que se adquiere senal de multiples direcciones 
utilizando los mismos codigos de correccion focal y diferentes valores N A (k) 
5 para diferentes direcciones con los que se adquiere el primer foco. En este 
caso, el valor de la direccion inicial debe restablecerse con anterioridad a la 
adquisicion. 

El dispositivo APD (60), Figura 4, realiza las funciones de apodizacion y 
apertura dinamica. La apodizacion, que consiste en multiplicar la senal recibida 

10 en cada canal por un numero real positivo A k comprendido entre 0 y 1, mejora 
las caracteristicas del haz conformado y el contraste en la imagen al reducir los 
lobulos laterales. La apertura dinamica permite mantener una resolution lateral 
constante operando con un numero F# constante, que requiere que el tamano 
de la apertura activa aumente progresivamente en el tiempo. Ademas, la 

15 apertura dinamica facilita la aplicacion de la tecnica de Correccion Focal 
Progresiva en las regiones mas proximas al transductor. 

La Figura 9 muestra un diagrama de bloques del dispositivo APD (60). Para la 
funcion de apodizacion utiliza un multiplicador MULT (61) que obtiene el 
producto de cada uno de los valores de la secuencia e k por el coeficiente A k de 

20 apodizacion asignado al canal k. En una realization preferente, este coeficiente 
es un numero entero sin signo de a bits. Si los valores e k se expresan con w 
bits, de los a+w bits que componen el producto se mantienen los w bits mas 
significativos con redondeo en funcion del valor del primer bit ignorado y el 
signo de e kl incluyendose la logica correspondiente en el propio dispositivo 

25 MULT (61) para obtener la salida h k . La funcion realizada es: 



K = 



A 



(47) 



Debe hacerse la observacion de que la funcion de apodizacion puede 
habilitarse o no con sencillos circuitos de control no mostrados en la figura 9. 
Basicamente, la salida h k es la entrada e k si la funcion de apodizacidn esta 
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inhabilitada, respondiendo a la Ecuacion (47) en otro caso, funcion que se 
realiza con un simple multiplexor gobernado por la senal de habilitacion de la 
apodizacion. 

Para la funcion de apertura dinamica, se dispone un contador binario CNZ (62) 
5 que se carga inicialmente con el valor MZ(k) o numero de muestras nulas 
expresado por la Ecuacion (25) para el canal k. Este contador se acciona con 
un reloj cuya frecuencia es la nominal de muestreo, cs. Su operacion se habilita 
con la activacion de la senal HFD=1 que inicia el proceso de focalizacion 
dinamica y es proporcionada por el dispositivo HFOC (80), Figuras 4 y 8. Por 
10 su parte, la senal HFD-1 actua sobre la senal S del flip-flop FF (64) poniendo la 
senal z-1 que selecciona la entrada nula (63) al multiplexor MUX (65). 

En este estado, se inicia la adquisicion y procesamiento de senales en todos 
los canales. En aquellos en los que el contador CNZ (62) no haya alcanzado la 
cuenta terminal (senal fc), las muestras de salida en f k seran de valor nulo, por 

15 lo que el sistema opera como si el elemento asociado no estuviera active Tras 
MZ ciclos del reloj de muestreo, durante los que se han producido otras tantas 
muestras nulas sobre f kl se alza la serial tc=1 que pone a cero el flip-flop FF 
(64) y habilita la salida de los valores h k obtenidos por el multiplicador (61). De 
este modo se proporciona un numero programable de muestras nulas al 

20 comienzo de cada adquisicion en cada canal, lo que permite programar las 
caracteristicas de la apertura dinamica. La secuencia f k es: 

f k = {0, 0, ...,0, A k e k (NZ+1), A k e k (NZ+2), ...} (45) 

donde e k (i) es la muestra / de la secuencia e k y que, por consiguiente, se 
corresponde con la muestra / de cualquier otra secuencia. De este modo se 
25 mantiene el orden para realizar posteriormente la composicion coherente de 
todas las senales, sean o no nulas. 

Otro de los aspectos importantes de la presente invencion es su capacidad de 
realizacion modular, cuyo concepto se muestra en la Figura 10. Una realizacion 
30 preferente agrupa en un modulo (50) todas las funciones de procesamiento de 
senal para un subconjunto de P S elementos del array, asi como otras 
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funciones de control global. En particular, un modulo incluye la memoria MJ 
(59) que proporciona los codigos de alejamiento focal J, junto con los circuitos 
necesarios para su lectura y carga inicial. Ademas, cada modulo (50) contiene 
S submodulos (40), donde cada submddulo procesa las senales de P 
elementos, y un combinador (55) que permite enlazar multiples modulos en 
cadena. 

Internamente, los S submodulos (40) que conforman un modulo, se encadenan 
para realizar la composicion coherente de los P S canales en cascada a traves 
de los FIFOs (52). Esta estructura facilita la operacion en paralelo de todos los 
canales de procesamiento digital, donde los FIFOs compensan 
automaticamente las diferencias en Ids instantes de adquisicion de serial en los 
diferentes canales. 

Para ello, cada uno de los elementos (11) del array (10) proporciona la serial 
recibida a uno de los P procesadores elementales, PE, (30) integrados en un 
submodulo (40). Cada procesador elemental PE realiza el acondicionamiento 
de senal y digitalizacion en el instante de llegada de la senal al elemento que le 
corresponde, siguiendo la metodologfa descrita anteriormente, e incluye los 
dispositivos que fueron mostrados con detalle en las Figuras 4 a 8. 

Cada PE proporciona a su salida una secuencia ordenada de valores 
correspondientes a las senales recibidas por el elemento asociado desde cada 
uno de los focos o muestras intercaladas. Estas salidas se suman (41) para 
realizar la composicion coherente de las P senales que corresponden a un 
submddulo. 

Por su parte, el dispositivo (42) suma los resultados obtenidos localmente en el 
submodulo con los que le haya colocado el precedente en su cola FIFO (52). 
De este modo, a la salida de cada submodulo se obtiene una secuencia 
ordenada que corresponde a las senales focalizadas de sus P elementos mas 
los de los submodulos que le preceden en la cadena, en el orden de llegada de 
senales desde cada foco. El primer submodulo de la cadena recibe un valor '0' 
(51) en sustitucion del FIFO (52), considerando que no existe un submodulo 
precedente. 
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A su vez, el combinador (55) suma en (58) los resultados obtenidos por este 
modulo y que se almacenan en el FIFO (56) con los que el modulo anterior ha 
colocado en la cola FIFO (57). Como en ambos casos los resultados son 
secuencias ordenadas de valores correspondientes a focos sucesivos, la 
5 secuencia de salida del sumador refleja este mismo orden. El encadenamiento 
de L modulos (50), cada uno de ellos con S submodulos (40) con P 
procesadores elementales en cada submodulo proporciona un sistema de 
focalizacion dinamica con la tecnica de Correccion Focal Distribuida para 
N=L S Pelementos. 

10 Es, por tanto, una estructura que presenta diversas ventajas: 

• Es modular, facilitando el desarrollo de circuitos y mayor flexibilidad de 
aplicacion. 

• Procesa en paralelo, lo que permite alcanzar elevadas tasas de calculo. 

• Tiene una capacidad ilimitada de expansion desde el punto de vista logico. 

15 La Figura 11 muestra una de las posibles realizaciones del sumador (41) 
integrado en un submodulo (40). Este sumador recibe las P secuencias e? a ep 
correspondientes a los valores de las muestras adquiridas por los P 
conversores A/D pertenecientes al sub-modulo. En el caso mostrado en la 
figura es P=8, siendo la estructura generalizable a cualquier otro valor, 

20 preferiblemente potencia de 2. 

Las muestras que corresponden a cada foco en las secuencias e 1 a e P se 
obtienen en forma ordenada, esto es, la primera muestra obtenida por los P 
conversores A/D corresponde al primer foco, la siguiente al segundo foco y asf 
sucesivamente. Este razonamiento es extensivo aplicado a las muestras 
25 obtenidas entre focos, con o sin la opcion de Correccion Focal Progresiva 
Variable. Sin embargo, los instantes de muestreo no son simultaneos en todos 
los canales, sino que es posible que un canal adquiera varias muestras 
sucesivas mientras que otros esperan a que las senales correspondientes 
lleguen al elemento asociado. 

30 Este desfase se corrige automaticamente con las memorias FIFOs distribuidas 
(43), (45) y (47) junto con los sumadores (44), (46) y (48). El FIFO (49) realiza 
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la misma funcion pero con referenda a las secuencias de muestras de dos 
submodulos. 

Los FIFOs se asocian en parejas, cuyas salidas se suman cuando ambos 
disponen de, al menos, una muestra almacenada. Esta situacion se detecta 
5 con un sencillo circuito de control que ejecuta una lectura simultanea en ambos 
FIFOs cuando no estan vacfos y escribe el resultado obtenido por el sumador 
correspondiente en el FIFO conectado a su salida. Cada lectura libera una 
posicion de memoria en los FIFOs, que puede ser ocupada con la escritura de 
una nueva muestra. Cada escritura ocupa una nueva posicion de memoria en 
10 el FIFO correspondiente. Las lecturas en los FIFOs asociados son simultaneas, 
pero no asf las escrituras que dependen de los instantes de muestreo en los 
respectivos conversores A/D, a su vez derivados del instante de llegada de la 
serial desde cada foco al elemento. 

La capacidad H de cada FIFO (43) debe ser suficiente para almacenar las 
15 muestras adquiridas hasta que se disponga de una muestra en el canal 
asociado, en cuyo momento opera el sumador (44), liberando la posicion 
ocupada por las muestras correspondientes y transfiriendo el resultado de la 
suma al FIFO (45). 

Analogamente, los FIFOs (45) deben tener una capacidad 2 H, los FIFOs (47) 
20 una capacidad 4 H, y asf sucesivamente. De este modo, la cantidad de 
memoria distribuida en forma de FIFOs requerida por un submodulo con P 
entradas es: 

HP + 2- HP /2 + 4-H*P/4 + ... = H-P-\og 2 P (49) 

En el caso en que se utilizara un unico FIFO por cada canal, su capacidad 
25 deberia ser de H (P-1) posiciones, con lo que el total de memoria requerida 
seria H P (P-1). Asf, por ejemplo, si H=16 y P=8, segun el metodo preferido 
cuya estructura se muestra en la Figura 6 se requieren 384 posiciones de 
memoria distribuida, mientras que en el caso de utilizar un unico FIFO por cada 
canal serian necesarias 896 posiciones. 

30 El FIFO (49), por su parte, almacena transitoriamente los resultados obtenidos 
por el submodulo actual hasta que el sumador (42) de la Figura 10 pueda 
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liberar una posicion cuando el FIFO (52) disponga de un dato correspondiente, 
entregando en tal caso el resultado al submodulo siguiente. 

Breve descripcion del contenido de las figuras. 

5 La Figura 1 muestra la geometria utilizada como ejemplo, con un array 
compuesto por N elementos numerados 1, 2, 3, ... k, AM y A/, sobre el eje x. 
Los focos se situan a intervalos AR a partir de una distancia minima R 0 . Se 
muestra el segmento entre un foco situado en coordenadas polares (R, Q) y el 
elemento k del array situado en coordenadas cartesianas (x k , 0), cuya longitud 
10 esL k (R,Q). 

La Figura 2 muestra la estructura general de un sistema de conformacion digital 
de haces en recepcion generico. Se compone de un array (10) con N 
elementos (11) y un canal de procesamiento de senal por cada elemento. Cada 
canal contiene circuitos de amplificacion y acondicionamiento de serial AAS 
15 (12) y un digitalizador A/D (13). Las salidas de los N digitalizadores son 
procesadas por el conformador CONF (14). 

La Figura 3 muestra la estructura utilizada frecuentemente para realizar la 
composicion coherente de sefiales segun metodos conocidos. Consta de un 
sistema de retardo grueso (22) y otro fino (23), cuya operacion se gestiona 
20 desde un controlador local (24) que utiliza el contenido de una memoria de 
focalizacion dinamica (25). Las N salidas retardadas f 1t f 2 , fa se suman en 
(20) para obtener la secuencia resultante r. 

La Figura 4 muestra la estructura de un canal de procesamiento digital 
preferentemente utilizado por el sistema objeto de esta patente. Contiene el 
25 elemento de amplificacion y acondicionamiento de senal AAS (12), el 
digitalizador A/D (13), una memoria FIFO (43) que almacena transitoriamente 
las muestras adquiridas, un generador de reloj de muestreo (70) que constituye 
un elemento fundamental de la presente invencion y un modulo de habilitacion 
de la focalizacion dinamica y control, HFOC (80). 

30 La Figura 5 muestra la estructura de principio para realizar el generador del 
reloj de muestreo (70) con codigos de correccion focal O expresados con 1 bit 
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en la memoria MEM (75). El contador CNT (76) establece la direcci6n de 
acceso a esta memoria. El reloj de muestreo CK k , propio del canal k, se genera 
mediante logica distribuida en el registro REG-A (71a), el registro de 
desplazamiento de longitud programable SHR-A (72), el flip-flop FFA (73), el 
5 multiplexor MUX (74) y el flip-flop FFB (77). 

La Figura 6 muestra la estructura de principio para realizar el generador del 
reloj de muestreo (70) con c6digos de correccion focal Q expresados con 
multiples bits en la memoria MEM (75). El contador CNT (76) establece la 
direccion de acceso a esta memoria. El reloj de muestreo CK k , propio del canal 
10 k, se genera mediante logica distribuida en el registro REG-A (71a), el registro 
de desplazamiento de longitud programable SHRA (72), el registro de 
desplazamiento de longitud programable SHRB (78) y el conjunto de inversores 
(75a). 

La Figura 7 muestra la estructura de principio para realizar el generador del 

15 reloj de muestreo (70) con c6digos de correccion focal Q expresados con uno o 
mas bits en la memoria MEM (75). El contador CNT (76) establece la direccion 
de acceso a esta memoria. En este caso puede insertarse un numero fijo o 
variable de muestras entre focos, segun el contenido del registro-contador 
RCM (71 b) controlando la variable ng que establece la longitud del registro de 

20 desplazamiento SHRB (78). La variable ng, de 1 bit, es proporcionada por el 
control CTRL (90) mostrado con mas detalle en la parte inferior de la figura. 
Esta unidad consta de un contador de modulo programable CMOD (91), cuyo 
contenido y el valor del codigo de correccion focal Q, forman la direccion a una 
memoria local MDM (92). La salida g de esta memoria es invertida por (93) 

25 para obtener la serial ng. En este caso, existe la opcion de proporcionar el 
codigo de alejamiento focal J desde la misma memoria MEM de codigos de 
correccion focal (75), como se muestra con trazos de puntos o, 
preferiblemente, desde el exterior y a traves de una memoria FIFO (79). Por lo 
demas, el reloj de muestreo CK k del canal k se produce mediante la logica 

30 compuesta por el registro REG-A (71a), el registro de desplazamiento SHRA 
(72) y el registro de desplazamiento SHRB (78). La unidad de control CTRL 



proporciona, ademas, la salida cef utilizada para habilitar el conteo de focos en 
la unidad HFOC (80) mostrada en las Figuras 4 y 8. 

La Figura 8 muestra esquematicamente el principio en que se basa la unidad 
HFOC (80) que realiza el control general de la adquisicion. Consta de un 
5 contador CNT-A (82) que se carga mediante Ida con el valor N A1 , se habilita 
con la serial cea y se acciona con el reloj de baja frecuencia cb. La cuenta 
terminal tea actua sobre la entrada S del flip-flop FF (83). Su salida aj se 
retrasa una cantidad de ciclos de reloj maestro c expr^sada por N A2 mediante 
el registro de desplazamiento SHRG (84). La puerta AND (81) produce un nivel 
10 alto en la salida HFD cuando sus dos entradas estan en alto. El flip-flop (83) se 
repone a 0 cuando se alza su entrada R al activarse la serial tcf. 

La senal tcf es la cuenta terminal del contador CNT-F (86) que puede 
inicializarse con un numero de focos N F determinado al activarse la senal Idf, se 
habilita con la senal cefy se acciona con el reloj de baja frecuencia cb. Por su 
15 parte, el contador CNT-S (88) cuenta las muestras adquiridas, activando su 
salida FIN cuando se alcanza el numero programado N s . Esta programacion se 
realiza activando la senal Ids, el contador se habilita con ces y se acciona con 
el reloj de baja frecuencia cb. 

La Figura 9 muestra esquematicamente el principio en que se basa la unidad 
20 APD (60) que realiza las funciones de apodizacion y apertura dinamica. Consta 
de un multiplicador (61) que obtiene el producto de la entrada e k por el del 
coeficiente A k y entrega su resultado h k a un multiplexor MUX (65). Asimismo, 
dispone de un contador CNZ accionado por un reloj cs y habilitado por la senal 
HFD cuyo contenido inicial es M z , esto es, el numero de muestras nulas a 
25 aportar por este canal para realizar la funcion de apertura dinamica. Al alzarse 
la senal HFD se pone en alto la senal z de salida del flip-flop FF (64), 
seleccionandose la salida '0* entregada por la unidad (63) por el multiplexor 
MUX (65), que la proporciona como salida f k . Cuando se alcanza la cuenta 
terminal tc se activa la entrada R del flip-flop FF (64), pasando a nivel bajo su 
30 salida z con lo que el multiplexor MUX (65) selecciona la salida h k del 
multiplicador MULT (61), entregandola por la salida f k . 
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La Figura 10 muestra esquematicamente la arquitectura modular del sistema 
objeto de esta patente. Un modulo (50) contiene multiples submodulos (40), 
cada uno de ellos con varios procesadores elementales PE (30) cuyas entradas 
son las senales recibidas por los elementos (11) del array (10). Las salidas de 
los procesadores elementales son sumadas en (41) y la salida de este sumador 
se utiliza como uno de los sumandos de una cadena de sumadores (42) y 
FIFOs (52). El primer submodulo sustituye este FIFO por una entrada nula (51). 
Los resultados proporcionados por el sumador (42) del ultimo submodulo se 
almacenan transitoriamente en un FIFO (56) para ser, a su vez, sumados en 
(58) con los resultados de un modulo precedente que han sido almacenados en 
el FIFO (57) en el sumador (58). Adicionalmente, el modulo contiene una 
memoria MJ (59) que proporciona los codigos de alejamiento focal J a los 
diferentes procesadores elementales PE (30). 

La Figura 11 muestra esquematicamente el principio en que se basa una 
realizacion preferida del sumador (41), con un ejemplo para 8 canales, donde 
multiples FIFOs (43) almacenan transitoriamente las muestras obtenidas en 
cada procesador elemental e k . Las salidas de los FIFOs (43) son, a su vez, 
sumadas por pares en los sumadores (44), almacenandose el resultado en los 
FIFOs (45). Las salidas de estos FIFOs son sumadas en (46) y el resultado se 
almacena en los FIFOs (47). Analogamente, las salidas de estos FIFOs son 
sumadas en (48), almacenandose los resultados finales en el FIFO (49), cuya 
salida r representa la suma de las muestras e 1 a e 8 entregadas por los 8 
procesadores elementales que constituyen un submodulo en este ejemplo. 

La Figura 12 muestra un esquema de un ejemplo realizado de un modulo (200) 
del sistema objeto de esta patente. Dispone de conectores (120a) y (120b) por 
los que se introducen las senales analogicas que son distribuidas por los buses 
(121a) y (121b) a los conversores A/D dobles (101). Las muestras 
proporcionadas por estos conversores son introducidas por (106) en una FPGA 

(102) que realiza las funciones de procesamiento digital correspondientes a un 
subm6dulo (100), entre las que destacan la generacion de los relojes de 
muestreo (105) y el suministro de resultados al siguiente submodulo por el bus 

(103) . El ultimo submodulo en la cadena suministra los resultados obtenidos a 
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una FPGA (111) por el bus (104) para ser combinados con los resultados 
procedentes de otro modulo por el bus (113), suministrandose los datos 
resultantes por los buses (112) o (114). Estos buses estan presentes en el 
conector (110). A su vez, la FPGA (111) entrega senales de control (115a) y 
5 (115b) a los circuitos de amplificacion y acondicionamiento analogico a traves 
de los conectores (130a) y (130b), respectivamente. 

La Figura 13 muestra las medidas experimentales realizadas sobre el modulo 
descrito en la Figura 12, mostrando el error absoluto en ns para cada uno de 
los 32 canales. 

10 EJEMPLO DE REALIZATION DE LA INVENCION: 

En todos los casos descritos se obtiene el reloj de muestreo de cada canal 
utilizando muy pocos recursos hardware, fundamentalmente registros, 
contadores y una memoria para almacenar los codigos de correccion focal, 
realizandose la adquisicion de las senales con una alta resolucion temporal. 

15 Debido a la alta eficiencia de la metodologla descrita en cuanto a utilizacion de 
recursos hardware y, a diferencia de otras realizaciones, las tecnicas CFP y 
CFPV permiten integrar multiples canales de procesamiento en un unico 
circuito VLSI, junto con otras funciones descritas, incluso en dispositivos de uso 
general como son las FPGAs. 

20 Asf, en una realizacion demostrativa de esta invention, un modulo contiene 4 
submodulos, cada uno de ellos con 8 procesadores elementales. Este pequeno 
sistema puede operar con una apertura activa de hasta 32 elementos en 
barrido lineal o sectorial con apertura dispersa y arrays de 128 elementos. Con 
4 modulos se configura un sistema que opera con aperturas densas o 

25 aperturas activas de hasta 128 elementos. El sistema es escalable, anadiendo 
32 elementos activos por cada nuevo modulo instalado. 

La Figura (12) muestra graficamente la configuration del modulo, realizado en 
una tarjeta de circuito impreso de 220x100 mm. En la figura se han omitido los 
circuitos correspondientes al preprocesamiento analogico de senales, con 
30 posible excitation de los transductores, asi como los de posprocesamiento que, 
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en esta realization particular, se basan en un computador integrado al que se 
conecta una Unidad de Control e Interfaz (UCI) con el sistema descrito. 

El modulo (200) recibe las senales analogicas de 32 canales proporcionadas 
por los circuitos de amplification y acondicionamiento por los conectores (120a) 
5 y (120b). Contiene cuatro submodulos (100) y un combinador realizado en la 
FPGA (111). El combinador recibe senales procesadas por un modulo anterior 
en la cadena por el bus (113) y envfa los resultados al que le sucede por el bus 
(114). Los parametros se programan desde la UCI a traves del bus de control 
(112) por el que, tambien, el ultimo modulo de la cadena entrega los resultados 
10 finales. Asimismo, la FPGA (111) entrega un conjunto de senales digitales por 
los conectores (130a) y (130b) que son utilizadas por la unidad de 
procesamiento analogico. Por lo demas, esta estructura refleja fielmente la 
arquitectura presentada en la Figura 10. 

Las senales analogicas introducidas por los conectores (120a) y (120b) se 
15 distribuyen a los submodulos mediante los buses (121a) y (121b), 
respectivamente, donde son digitalizadas por los conversores A/D (101). La 
salidas digitales de los conversores A/D se introducen en las FPGAs (102) 
mediante los buses (106). Los relojes de muestreo de dichos conversores (105) 
se generan en las FPGAs (102) con la tecnica de Correction Focal Progresiva 
20 descrita en esta memoria. 

Por su parte, los submodulos se conectan entre si mediante los buses 
bidireccionales (103). Junto al bus (104) que conecta el ultimo submodulo al 
combinador (111), estos buses se utilizan en un sentido para la programacion 
de parametros desde el de control (112) y en el sentido opuesto para 
25 suministrar resultados a los submodulos posteriores en la cadena o al 
combinador (111). 

El combinador realiza la suma coherente de los resultados obtenidos en el 
modulo actual a traves del bus (104) con los que recibe de un modulo anterior 
por el bus (113), suministrando los resultados de este proceso por el bus (114) 
30 al modulo siguiente o por el bus (112) a la Unidad de Control e Interfaz. Una u 
otra option son programables. 
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En esta realizacion se han utilizado 16 conversores dobles de 10 bits, tipo 
AD9218 de Analog Devices para los dispositivos (101), 4 FPGAs tipo XC2S200 
de Xilinx Inc. para los dispositivos (102) y una XC2S150 de Xilinx Inc. para el 
dispositivo (111). 

5 El reloj del sistema es de 40 MHz, velocidad a la que operan sincronamente las 
comunicaciones entre submodulos, modules y con la UCI. Internamente, parte 
de los circuitos tambien operan a esta frecuencia o una fraccion que determina 
la de muestreo. Sin embargo, en las FPGAs (102) se multiplica por 4 la 
frecuencia del reloj del sistema para obtener el reloj maestro de 160 MHz. La 
10 resolucion temporal es, pues, T x = 6.25 ns. y el error de muestreo esta 
teoricamente acotado a I e I < 3.125 hs. Con esto, pueden procesarse senales 
analogicas de hasta 10 MHz con [i > 16 y hasta 20 MHz con \i > 8, siendo \x la 
relacion entre el periodo de la serial y la resolucion temporal que afecta al 
rango dinamico y contraste de las imagenes. 

15 Cada FPGA (102) contiene los circuitos necesarios para realizar todo el 
procesamiento digital de 8 canales, incluyendo: 

• Apodizacion, con coeficientes 0 < A k < 255/256 o deshabilitada (A k = 7, V /c). 

• Apertura dinamica, con 0 < MZ k < 2 12 . 

• Deflexion del haz con foco fijo, con 0 < NA k < 2 14 . 

20 • Focalizacion dinamica segun la tecnica de Correccion Focal Progresiva 
expuesta, con 4096 codigos de correccion focal de 1 bit por cada canal. 

• Programacion del numero de muestras a adquirir, con 2<N S < 2 17 . 

• Programacion del numero de muestras entre focos, con 4 < m < 64. 

• Programacion del numero de focos por adquisicion, con 0 < N F < 2 12 . 

25 • Programacion de la frecuencia de muestreo de salida, con 2.5 < f s < 40 
MHz. 

• Circuitos de gestion de ahorro energetico y habilitacion de funciones. 

• Circuitos de programacion de parametros y procesamiento de resultados. 
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• Circuitos para la verification del modulo, integrando un generador digital de 
serial. 

Por su parte, el combinador contenido en la FPGA (111) tiene circuitos para 
5 realizar: 

• La gestion del modulo, transfiriendo valores de parametros a los 
submodulos por (104) y (103). 

• La recopilacion de resultados procedentes del modulo anterior por (113). 

• La composition de estos resultados con los del propio modulo, recibidos por 
10 (104). 

• El suministro de resultados combinados al modulo siguiente por (114). 

• El suministro de resultados finales a la UCI a traves del bus (112). 

• Circuitos para compensar el offset de las senales. 

• Circuitos para controlar la unidad de pre-procesamiento analogico. 

15 

La Figura (13) muestra medidas experimentales relativas al error absoluto en 
ns. (ordenadas) del instante de muestreo en cada uno de los 32 canales (en 
abcisas) para un foco situado a una distancia de 50 mm., con un angulo de 
deflexion de 15°, velocidad de propagation c= 5900 m s' 1 y array lineal de 5 
20 MHz con elementos uniformemente distribuidos a intervalos d=0.8 mm. 

El maximo error teorico es de 3.725 ns., que se verifica en 28 de los 32 
canales. El mayor error en los 4 canales que superan el valor teorico se debe a 
diferencias en el trazado de pistas en el circuito impreso y en las FPGAs de 
este modelo experimental y, en cualquier caso, no es significativo. EI error 
25 cuadratico medio, mas determinante de la calidad de la imagen, es de 7.9 ns. 

Es evidente que pueden introducirse mejoras a la realization descrita puesto 
que existen dispositivos con mayor capacidad de integration. Asimismo, 
pueden realizarse otras modificaciones que afectan al tipo de conversion A/D, 
distribution de dispositivos, etc., ya que el modelo descrito no es mas que un 
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ejemplo demostrativo, sin que tales modificaciones supongan un cambio 
sustancial a los fundamentos de la metodologia y sistema que han sido 
descritos en la presente memoria. 



5 Para la aplicacion del metodo, considerese un caso de imagen ultrasonica 
sectorial en pulso-eco, con un transductor lineal de N elementos centrado en el 
origen, con un tamano de apertura D. Se pretende cubrir con la imagen un 
sector de -45° < 9 < 45°. 

Si f R es la frecuencia de la senal recibida, para verificar con suficiente margen 
10 el criterio de Nyquist, se elige una frecuencia de muestreo nominal fs = 4 f R . 
Ademas, para mantener un elevado rango dinamico y contraste en la imagen, 
se elige = f x /fR =16, siendo f x la frecuencia del reloj maestro. 

Considerando el trayecto en ida y vuelta del pulso ultrasonico, el intervalo entre 
muestras es: 

15 AR = cT s /2 (50) 

Sustituyendo el valor AR dado por la Ecuacion 1 resulta, con m = 1: 



(51) 



T 
1 x 

Cuando se insertan m muestras entre focos, el valor de v = 4m. Aplicando 
ahora la Ecuacion 18, con \x k \ < D/2, donde D es el tamano de la apertura, la 
20 distancia minima para aplicar la Tecnica de Correccion Focal Progresiva con 
apertura completa es: 



D 2 



mv -2a 

COS0 



2ylamv - a 2 



+ -~|sen0| (52) 



donde a=2 b -1, siendo b el numero de bits con los que se expresan los codigos 
de correccion focal, con la condicion mv > 2a. 

25 La Tabla II muestra los resultados de la Ecuacion (52) para b=1, 2 y 3 en 
funcion del numero de muestras (fijo) entre focos. 
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Debe observarse que el valor R</D representa el minimo valor del numero F# 
con el que se puede operar utilizando la funcion de apertura dinamica. En 
general, se acepta como valido un valor F# min < 1, que se obtiene con la 
mayoria de las combinaciones. 

Tabla II 

Relacion F# min = Rc/D para aplicar la TCFP con apertura completa (-45°< 9 < 
45°) 



m 


b = 1 


b = 2 


b = 3 


1 


0.6 






2 


0.8 


0.4 




4 


1 


0.6 


0.4 


8 


1.4 


0.9 


0.6 
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REIVINDICACIONES 

Esta patente descubre el metodo y sistema para conformacion de haces en 
recepcion con una nueva tecnica denominada de Correccion Focal Progresiva, 
con la opcion Variable. Segun dichos metodo y sistema, se realiza la 
5 composicion coherente de senales procedentes de fuentes situadas a lo largo 
de una direccion seleccionada, que llegan a un conjunto de N transductores 
receptores o elementos distribuidos espacialmente, formando un array. Un 
conformador de haces en recepcion realiza una composicion coherente de las 
senales recibidas que proceden de fuentes puntuales situadas a lo largo de 

10 dicha direccion. La composicion coherente de las senales originadas en cada 
punto y recibidas por cada uno de los N elementos del array proporciona una 
unica senal de salida que representa a la que se obtendria con un transductor 
receptor unico, del mismo tamano que el array, focalizado individualmente en 
todos y cada uno de los focos, tecnica que se conoce en el ambito como 

1 5 focalizacion dinamica con deflexion. 

Las senales recibidas son amplificadas y acondicionadas analogicamente antes 
de ser digitalizadas. La tecnica de correccion focal progresiva, con o sin la 
opcion variable, y el sistema que la realiza descritos en esta patente, se aplican 
a las senales digitalizadas por un conversor analogico-digital. Se denomina 
20 canal al conjunto de dispositivos que procesan la senal recibida por un 
elemento del array hasta su combinacion con las correspondientes a otros 
elementos. 

A diferencia de otros metodos que forman parte del estado del arte, el metodo y 
sistema objeto de esta patente no se basa en retrasar las senales recibidas por 

25 cada uno de los elementos del array, ni utiliza procesos de interpolacion. Antes 
bien, se fundamenta en muestrear las senales recibidas en los instantes en que 
llegan a cada uno de los elementos del array desde cada uno de los focos. El 
metodo garantiza que la composicion de muestras es coherente, esto es, solo 
se combinan las que corresponden al mismo foco y, ademas, no redundante, 

30 adquiriendo unicamente las que van a utilizarse en el proceso de composicion. 
Requiere un reloj de muestreo de las senales no uniforme y diferente para cada 
uno de los canales. 
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Por otra parte y, a diferencia de otras tecnicas en las que cada reloj de 
muestreo queda establecido por una secuencia de '1' y '0' en una memoria 
individual para cada canal o se calcula en tiempo real mediante circuitos 
digitales especfficos a partir de los valores de una diversidad de parametros 
5 programados durante la adquisicion, con el metodo y sistema objeto de esta 
patente los instantes de muestreo se determinan interpretando el contenido de 
tablas asociadas a cada canal, que se calculan y almacenan en memorias 
antes del inicio de la adquisicion de serial, haciendo innecesaria la carga de 
multiples parametros en tiempo de adquisicion. Ademas, la cantidad de 

10 almacenamiento requerida es limitada al indicar, en forma codificada, la 
correccion focal que debe aplicarse al reloj de muestreo para obtener el 
siguiente foco. Mas aun, cada codigo de correccion focal puede ser compartido 
por multiples muestras consecutivas. Finalmente, la arquitectura de la tecnica 
aquf patentada presenta una elevada regularidad, facilitando su integracion en 

1 5 dispositivos VLSI y su realizacion modular. 

Se reclaman, como propias, las siguientes reivindicaciones: 

1. Un metodo de conformation de haces en reception, denominado Tecnica 
de Correccion Focal Progresiva, para aplicaciones en las que una 
diversidad de N transductores, agrupados en un array, reciben energfa 

20 vibratoria desde focos F h i=0, 1, 2,..., situados a lo largo de una direction 
programable, caracterizado porque muestrea la senal recibida en cada 
canal en el instante de llegada desde cada foco con un error absoluto 
inferior a medio periodo de un reloj maestro de periodo T x , de forma que las 
muestras consecutivas e kh /=0, 1, 2, ... obtenidas en cada canal k 

25 correspondan a los valores consecutivos de la senal originada en cada foco 
F u i=0, 1, 2,..., con lo que cada muestra de la serial resultante r h i=0, 1, 2, ... 
se realiza mediante la suma de las N muestras e kh k=1, 2, .., N : 

N 

r i=S c « i = 0,l,2,... 
o de una combinacion lineal como: 



N 



i = 0, 1, 2 



9 • •• 



donde los coeficientes A k son valores fijos o programables, operacion que 
se denomina suma o composicion coherente de las N senates, 

Un metodo de conformacion de haces en recepcion segun reivindicacion 1 
caracterizado, ademas, porque para el canal receptor correspondiente al 
elemento k del array, que incluye circuitos para el acondicionamiento y la 
amplification de senal, el instante de muestreo que corresponde al primer 
foco F 0 se determina, en tiempo no real, a partir del numero de ciclos N A (k) 
del reloj maestro que median entre un origen de tiempos elegido con la 
restriction de que todos los tiempos sean positivos y el instante de llegada 
de la serial del foco F 0 al canal /c, de modo que si el tiempo de vuelo de la 
senal correspondiente al foco F, al elemento k es T,(k), calculado a partir de 
la geometria del sistema y de la velocidad de propagation de las ondas en 
el medio, el valor N A (k) viene dado por: 



donde [ ]n representa la funcion de redondeo al entero mas proximo, que 
determina N A (k) con un error absoluto inferior a T x /2, funcion que puede 
realizarse mediante un contador de N A (k) ciclos de reloj maestro, 
generandose un pulso de adquisicion de senal en el canal ken el momenta 
en que dicho contador rebosa; tal y como se ha descrito en la memoria, 
existen diversas formas de realizacion fisica, prefiriendose una en la que el 
contador se desdobla en dos partes, una de pequeno tamano accionada 
por el reloj maestro y la principal por un reloj cuya frecuencia es un 
submultiplo de la del reloj maestro 

Un metodo de conformacion de haces en recepcion segun reivindicaciones 
anteriores caracterizado, ademas, porque para otros focos F h i > 0, el 
instante de muestreo en el canal k se determina en tiempo real a partir del 
valor Q k f de b de bits, almacenado en una memoria de correccion focal 



T 0 (k) 



T 

x 
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asociada a dicho canal k que codifica el adelanto a aplicar al valor nominal v 
de periodos del reloj maestro que hay entre dos focos consecutivos, donde 
el codigo Q ki se calcula en tiempo no real a partir de la diferencia AT ki entre 
el instante de llegada al elemento k de las senales procedentes del foco i-1 
y el foco / como: 



ft,-=v- 



n 



donde [ ]u representa la funcion de redondeo al entero mas proximo yves 
el intervalo nominal entre focos, expresado en periodos del reloj maestro, 
programandose consecutivamente los codigos Q ki que corresponden a 
focos sucesivos F h i=1, 2, 3, ... en una memoria de correccion focal 
asociada al canal k cuyo contenido se lee en tiempo de adquisicion de 
serial para producir el adelanto variable en el instante de muestreo para 
cada foco F, en funcion del valor de Q ki mediante circuitos digitales de los 
que en la memoria descriptiva se han dado varios ejemplos, basados en 
registros de desplazamiento de longitud programable, siendo posibles otras 
alternativas, como el uso de contadores programables, sin que por ello se 
modifique la idea base expuesta en esta reivindicacion. 

Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores caractertzado porque permite insertar un numero 
fijo o programable m de muestras entre focos, por medio de un contador 
programable modulo N M -m que se incrementa cada vez que se adquiere 
una nueva muestra y, cuando retorna a cero, se incrementa la direccion de 
acceso a la memoria de correccion focal, facilitando que un total de m 
muestras compartan el mismo codigo de correccion focal Q ki , lo que permite 
aumentar ia eficiencia de la codificacion que puede ser de una fraccion de 
bit por muestra adquirida, manteniendo el error del instante de muestreo 
basicamente acotado dentro de medio periodo del reloj maestro al 
disponerse de circuitos digitales que distribuyen uniformemente el adelanto 
Q k i entre las m muestras que se adquieren entre focos, circuitos que se 
basan en una tabla en la que se programan, antes de la adquisicion y para 



cada posible valor Q ki < 2 b -1, un total de m bits g que indican, para cada una 
de las m muestras, si debe adelantarse (g=7) o no (g=0) en 1 ciclo de reloj 
maestro el instante de adquisicion de dicha muestra, constituyendo un caso 
particular de la metodologla expuesta en las reivindicaciones anteriores. 

5 5. Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores caracterizado porque puede variarse la distancia 
entre focos mediante la adquisicion de un numero variable de muestras m 
entre focos que se incrementa a partir de un valor inicial m 0 cuando el valor 
de un codigo de alejamiento focal, J h de 1 bit, asociado al foco F/ es J t - 1, 

10 opcion denominada tecnica de Correccion Focal Progresiva Variable, que 
aumenta la eficiencia de uso de la memoria de correccion focal al aumentar 
con el tiempo la distancia entre focos, siendo el valor de J, comun a todos 
los canales para cada valor de /, pudiendo tabularse de forma totalmente 
centralizada para todos los canales, parcialmente centralizada afectando a 

15 un conjunto de canales, como se realiza preferentemente, o individualmente 
en cada canal, lo que permite establecer un balance entre la flexibilidad de 
aplicacion y la eficiencia de memoria; para su realizacion basta sustituir el 
registro del valor constante m asociado a cada canal k por un contador que 
se carga con el valor inicial m 0 antes de iniciar cada adquisicion y, en 

20 tiempo de adquisicion, dicho contador se incrementa cada vez que se 
reciben senales desde un nuevo foco F, y el codigo Jj=1 . 

6. Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores caracterizado porque dispone de contadores 
para determinar el numero de focos N F y de muestras N s a adquirir, 

25 teniendo en cuenta que cada foco puede incluir la obtencion de m muestras, 
siendo m constante o variable, de forma que cuando el contador N F rebosa, 
se bloquean los mecanismos de la correccion focal progresiva y de 
actualizacion de los codigos de alejamiento focal, congelandose, ademas, la 
direccion de memoria que accede a los codigos de correccion focal en cada 

30 canal, aunque la adquisicion de senales puede continuar hasta que el 
contador de muestras a adquirir N s alcance el valor terminal, facilitando que 
una unica memoria por cada canal pueda almacenar los codigos de 



correccion focal para multiples angulos de deflexion, multiples 
configuraciones de la apertura activa o una combinacion de ambas, 
pudiendo programarse en cualquier caso la direccion inicial de la memoria 
de correccion focal antes de cada adquisicion u obtenerse de forma 
5 automatica si la secuencia de adquisiciones se corresponde con el orden en 
que se han almacenado los codigos de correccion focal para la actual 
aplicacion. 

7. Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores caracterizado porque opcionalmente puede 

10 operar con apertura dinamica, activando un contador asociado a cada canal 
k, en el que se programa un numero N z (k) de muestras de salida de valor 
nulo que se producen a partir del comienzo de la adquisicion para su 
combinacion con las muestras correspondientes de los restantes canales, 
de modo que, para estas N z (k) muestras, el resultado equivale a eliminar el 

15 elemento k en la composicion de las correspondientes muestras de salida, 
realizando efectivamente la funcion de apertura dinamica. 

8. Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores caracterizado, ademas, porque puede insertar 
una funcion de apodizacion o multiplicacion de los valores de cada muestra 

20 adquirida por una constante programable para cada canal. 

9. Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores caracterizado porque realiza la composicion 
coherente mediante una combinacion de colas o memorias FIFO 
distribuidas y sumadores, preferentemente formando una estructura en 

25 arbol en la que cada dos FIFOs se asocian a un sumador de 2 entradas. 

10. Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores caracterizado porque facilita su realizacion 
modular con una arquitectura flexible que permite la composicion coherente 
de un numero indefinido de senales desde el punto de vista logico, para lo 

30 cual, cierto numero de canales se agrupan para formar un sub-modulo, un 
numero de sub-modulos configuran un modulo y un sistema se configura 
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con cierto numero de modulos, donde todos los canales son iguales y 
realizan el procesamiento en paralelo de la informacion adquirida. 

11. Un metodo y sistema de conformacion de haces en recepcion segun 
reivindicaciones anteriores, que utiliza la metodologla y modelos 
constructivos expuestos en la memoria y figuras adjuntas. 
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